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Abstract 

Cli m ate simu la tions by 1 6 atmospheric general circula tion models (AGCMs) are compared on an aq ua-pla net, a 
wa ter-covered Earth with prescribed sea smface temperatme vm-y ing only in latitude. The idealised configura tion is 
designed to ex pose differences in the ci rculation simu lated by different models . Basic featu res of the aqua-planet 
climate are characterised by comparison with Emth. 

The models d isplay a w ide range of behaviour . The balanced com ponent of the tropospheric mean flow, and 
mid-latitude eddy covariances subject to budget constraints , vary relatively little among the models . In contrast , 
differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical unceltainty 
in modelled lower stratospheric temperatures persists in APE. 

Aspects of the circulation generated more directly by interactions between the resolved fluid dymunics and param- 
eterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter- tropical 
convergence zone (ITCZ) at the equator, with large variations in mean precipitation .The equatorial wave spectrum 
shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propa- 
gation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive 
than the equatorial Rossby modes, dominates in a few models or occm-s within an eastward propagating envelope 
in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies. 

The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top 
of atmosphere globa 1 net flux , dominated by differences in shmtwave reflection by clouds . A num ber of newly 
developed models , not optimised for Earth climate, contribu te to this. Possible reasons for differences in the 
optimised models are discussed. 

The aqua- planet configuration is intended as one com ponent of an experimental hiermchy used to evalua te 
AGCMs . Ill is comparison does suggest that the range of model behaviour could be better understood and reduced 
in conju nction with Earth climate simul ations . Controlled experi menta tion is req uired to explore indi vidual model 
behaviour and investigate convergence of the aqua-planet climate with increasing resolution. 

Keywords comparison of atmospheric general circulation models (GCMs); idealized model configuration; global 
energy buget; tropical wave spectnun; precipitation 


I. Introduction 

The Aqua-Planet Experimen t (APE) is a coordi- 
nated comparison of Atmospheric General Circula- 
tion Model (AGCM) simulations on a water-covered 
Earth . It compares idealized climates simulated by 
global AGCMs being developed and used for numer- 
ical weather prediction and climate research . The 
experimen t involves complete AGCMs but on an 
idealized planet with simplified lower boundary inter- 
actions, namely an ocean surface everywhere with no 
land, orography or sea ice. The sea surface tempera- 
ture (SST) is specified everywhere with a number of 
idealized analytic distributions designed by eale and 
Hoskins (2000a). 

APE was conceived by Neale and Hoskins as one 
component of a modelling hierarchy of increasing 
complexity, both of the models themselves and of the 
experimen tal configurations to which they are applied. 
The hierarchy has two distinct roles: an evaluation rol e 
in the development and testing of complex atmospheric 
models, and a conceptual role in 1 inking complex 
models with themy and observation. This context of 


APE is discussed in more detail by Blackburn and 
Hoskins (2013), so only a summary is given here. 

In the conceptual context, theory and more ideal- 
ised models provide constraints on the character of the 
global circulation expected in APE. The zonally aver- 
aged SST distributions in APE are broadly similar to 
Earth, so the equator-to-pole temperature difference is 
expected to give rise to a jet stream and storm track in 
mid-latitudes . In the tropics, the analytic model of Held 
and Hou ( 1980) predicts that a Hadley circulation with 
equatorial ascent must exist when the latitudinal profile 
of temp erature in equilibrium with the underlying SSTs 
is steeper than quartic, in order to maintain a therm al 
wind balance consistent with angular momentum 
conservation . The tropical SST profiles in APE were 
designed to include this limiting case . The CO TROL 
SST case considered in this paper is quadratic, so the 
tropical circulation is expected to consist of a Hadley 
circulation in each hemisphere with equatorial ascent. 

Previous modell ing studies using aqua-planet 
configurations have confirmed that the tropical circula- 
tion and the location of the inter-tropical convergence 
zone (ITCZ) produced by AGCMs do depend on SST. 
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The most systematic study, by Hess et al. (1 993), used 
a number of widely varyi ng tropical SST profi les and 
two different convective parameterizations i n a single 
model. Hess et al. and other studies found that the trop- 
ical circulation and ITCZ location also differ between 
models, with some models producing a double-ITCZ 
straddl ing the equator even when SST is peaked on 
the equator. Such a double-ITCZ was found in the 
fi rst modem aqua-planet simulation, by Hayashi and 
Sumi ( 1986). The modelled tropical circul ation can 
be affected not only by major modell ing choices, such 
as the type of convective parameterization (e.g., Hess 
et al. 1993), but also by more detailed choi ces such as 
model time -step (Will iamson and Olson 2003). ITCZ 
location can also be affected by disabl ing specific 
physical processes and interactions in models, such 
as the wind speed dependence of surface evaporation 
(Numaguti 1993; Chao and Chen 2004). 

These studies, in particular that of Hess et al., 
provided motivation for Neale and Hoskins (2000a) 
to propose a benchmark test suite of AGCM aqua- 
planet experiments. In this eval uation context, APE is a 
bridge in the modelling hierarchy between experiments 
with models ofreduced complexity that are used in the 
development of individual model components, and 
realistic simulations of Earth's climate using complete 
AGCMs, coordinated through the Atmospheric Model 
Intercomparison Project (AMIP, described in its orig- 
inal form by Gates 1992). 

Neale and Hoskins (2000b) used a single AGCM to 
characterise the climate for each of the APE SST distri- 
butions . The aims of the APE coordinated compar- 
ison described here are to provide a benchmark of 
current model behaviours and to stimulate research to 
understand the cause of inter-model differences. For 
this reason APE is sponsored by the World Meteoro- 
logical Organization (WMO) joint Commission on 
Atmospheric Science (CAS), World Climate Research 
Program (WCRP) Working Group on Numerical 
Experimentation (WGNE). 

In this and a companion paper (Williamson et al. 
20 I 3) we provide an overview of a variety of aspects 
of the aqua-planet simulations . Detailed analysis of 
specific aspects is presented in other papers in this 
Special Issue or left for future individual studies. This 
paper considers simulations based on the CONTROL 
SST distribution (defined below). The simulated 
responses to varying the latitudinal profile of SST are 
discussed in the companion paper, using experiments 
based on the other zonal ly symmetric SST distributions 
defined by Neale and Hoskins. More comprehensive 
diagnostics are presented in Williamson et al. (2012), 


a collection of model simulation statistics that wil 1 be 
referred to in the fol lowi ng as the APE ATLAS . 

2. APE experimental strategy 

2. 1 Sea suface temperature 

The A PE desi gn originated with Neale and 
Hoski ns (2000a) who specified experiments based on 
eight different SST fields, of which fi ve are zonally 
symmetric varying in latitude onl y and three impose 
tropi cal longitudinal perturbations on one of those 
zonally symmetric distri butions. The phi losophy of the 
design is described in Neale and Hoskins (2000a) and 
in Blackburn and Hoskins (2013) and is not repeated 
here. Only one of those 8 profiles is considered here 
which is zonally symmetric and labelled CONTROL. 
The specified zonally symmetric SST in °C is given by 
27 [1 -sin ^ (301/2)] for latitude ji between ±n/3, 
and 0 for I 0/ 1 2:: n 13 . Compared to present day Earth, 
this SST profile is more sharply peaked about the 
equator and has smaller equator-to -pole temperature 
difference because of the 0°C limit beyond 60° lati- 
tude. 

2 . 2 Detailed specification 

Certain model specifications were included in the 
APE design to obtain uniformity across the models. 
They are divided into Requirements and Recommen- 
dations. Complete details can be found in the APE 
ATLAS and on the APE website. [Available at http:// 
climate.ncas.ac. uk/ape/.] 

The insolation is perpetual equinoctial , symmetric 
about the equator, but including the diurnal cycle . The 
recommendation to achieve this is to modify the Earth 
orbit parameters , setting eccentricity and obliquity to 
zero, giving a circular equinoctial orbit. The distribu- 
tion of insolation is then independent of the calendar. 
The solar constant is 1365 W m- 2 . Ozone is taken from 
the annual mean climatology used in AMIP 11 (Wang 
et al. 1995; Liang and Wang 1996), symmetrized about 
the equator. C02 is to be set at 348 ppmv, as in AMIP 
II. Recommendations for the other well-mixed radia- 
tively active gases also follow AMIP II: CH 4 at 1650 
ppbv and N10 at 306 ppbv. Halocarbon concentrations 
should yield 0.24 W m- 2 radiative forcing. The use of 
an "equivalent " CO 2 is discouraged. It is also recom- 
mended to exclude radiatively active aerosols. Any 
aerosol specification for cloud condensation should use 
an oceanic distribution which is fixed in time, zonally 
symmetric and symmetric about the equator. Recom- 
mended values are provided for geophysical constants 
and parameters. 

It is recommended that the initial dry mass of the 
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Table I . Requested and derived diagnostic data 



DESCRIJ>TJON 

TIME AVERAGING 

DOMAIN 

GT 

Global Time-Series 

Daily 

Global, Area Average 

TR 

Transients 

6-hourly (year 3) 

Global/tropical lat-Jon 

SH 

Single-Level 2-D Means Multiple- 

3-yem-Jy 

Global lat-Ion 

ML 

Level 3-D Mean 

3 -yearly 

Global lat-Ion, P 17 1 

MF 

Multiple-Level Fluxes 

3 -yearly 

Zonal average, P 1 7 

PF2 

Parametrization Forcing 

3 -yearly 

Global lat-Ion, model- level 

TE 3 

Transformed Eulerian Means 

3-yem-Jy 

Zonal average, P I 7 

VB3 

Veltically Integrated Budgets 

3-yem-Jy 

Zonal average 


^17 denotes 17 WMO pressure levels, 10 1000 hPa. 
^Optional. 

3 Derived. 


atmosphere be equivalent to a global mean surface 
pressure of 10 1080 Pa. This is 101325 Pa minus 245 
Pa, which corresponds to a global moisture content of 
25 .006 kg m- 2 . It is desirable that dry mass be conserved 
throughout the integration , but this, of course, depends 
on the characteristics of the model dynamical core. 

2.3 Simulation period 

The simulations are run for 3.5 years, with the 
anal ysis made over the last 3 years, omitting the first 
6 months as spin-up. A model-simulated state, from 
either a previous aqua-planet integration or an earth- 
like simulation, should be used for the initial condi- 
tions. In independent tests of the spin-up in several 
models , aqua-planet climate equilibrated in a matter 
of weeks from these types of initial condition, so a 6 
month spin-up is considered adequate. Nevertheless, 
the modelling groups were instructed to check that 
equilibration was achieved during the discarded initial 
6-month period. Daily time-series of global averages 
were requested for a number of variables , which allows 
a gross check that equilibrium was reached . 

The book-keeping calendar can be a 365 or 360 day 
year, with natural or 30-day months respectively. Inso- 
lation does not follow the calendar, ote that the 3.5 
year integration length means that a realistic calendar 
can be used if integrations are started in March of a 
leap year. 

2. 4 Diagnostic data 

The APE proposal requested the submission to a data 
archive of a variety of data from the final three years of 
the simulations, in the categories listed in Table I. The 
diagnostic list is greatly expanded beyond the variables 
suggested by Neale and Hoskins (2000a). A complete 
list of variables in each category is included in the 


ATLAS and is available from the APE website. 

The required diagnostics comprise the top 5 catego- 
ries in Table 1. Global energy and moisture balances 
are analysed using the time series of global ly averaged 
single level variables (GT). These also allow simula- 
tions to be checked for residual drift. Temporal vari- 
ability of selected dynamical variables and tropical 
cloud and preci pitation are analysed using 6-hourl y 
time-series of global/tropical distributions of selected 
single-level variables (TR). The time average state is 
obtained from the 3-D atmospheric variables (ML) 
and 2-D single level variables (SH). The balance of 
processes that maintains the zonal average atmo- 
spheric state is obtained from the dynamical variances 
and covariances ("fluxes", MF) and the requested 
(optional) parameterization forcing (PF) diagnostics. 

The APE proposal also included specification of 
derived diagnostics that can be computed from the 
required and requested diagnostics . These comprise 
individual variables in each of the preceding catego- 
ries, plus two additional categories of time average 
data : zonal averages of Transformed Eulerian Means 
(TE) and vertically integrated budgets (YB). 

The standard diagnostics were considered necessary 
for the adequate analysis and comparison of the exper- 
iments. The list is an amalgamation of requested diag- 
nostics taken from the AMIP II standard model output, 
the WGNE standard diagnostics of mean climate 
and selected diagnostics that have proved useful in 
analyzing previous aqua-planet and dynamical core 
experiments. Many, but not all , of those diagnostics 
along with additional ones that have proven informa- 
tive are included in the ATLAS. A selection of those 
is discussed here to illustrate the general characteris- 
tics of the aqua-planet cl imate and the variation among 
models. 
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Table 2. Participating models 



GROUP SYMBOL 

LOCATION 

MODEL 

OPTIMISED TOA? 1 

1 

AGU 

Japan (consortium) 

AFES 

No 

2 

CGAM 

Reading, UK 

HadAM3 

Yes 

3 

CSIRO 

Aspendale , Australia 

CCAM 

No 

4 

DWD 

Offenbach /Mainz, Germany 

GME 

Yes (NWP) 

5 

ECM-CY29 

Reading, UK 

IFS cy29r2 

No 

6 

ECM-CY32 

Reading, UK 

IFS cy32r3 

No 

7 

FRCGC 

Y okohama , Japm 1 

NICAM 

No 

8 

GFDL 

Princeton, USA 

AM2.1 

Yes 

9 

GSFC 

Maryland, USA 

NSIPP-1 

No 

10 

K1JAPAN 

Japan (collaboration) 

CCSR/N1ES 5.7 

Yes 

11 

LASG 

Beijing, China 

SAMIL 

No 

12 

MIT 

Cambridge, USA 

MIT-GCM 

No 

13 

MRI 

Tokyo, Japan 

MR1/JMA98 

No 

14 

NCAR 

Boulder, USA 

CCSM-CAM3 

Yes 

15 

UKMO (N48) 

Exeter, UK 

pre-HadGAM 1 

Weakly 2 

16 

UKMO (N96) 

Exeter, UK 

pre-HadGAM 1 

Weakly 2 


1 Was the top of atmosphere radiative balance optimised for present day Earth climate? 

2 During the development phase ofHadGEMl the TOA fluxes of AM1P nms were monitored to check that 
they did not widely diverge from balance but were not actively tuned. 


All standard diagnostic data received from the APE 
modell ing groups are available from a Data Archive at 
the APE website . 

2. 5 Data processing 

Diagnostic data were supplied on a regular lati- 
tude-longitude grid unique to each model. In most 
cases this is the model’s native grid, in the case of 
spectral transform models the transfonn Gaussian 
grid. Data from models with non-regular grids were 
first interpolated to a regular grid by each modelling 
group . In addition, ECMWF interpolated from the 
model’s transform grid to a uniform 2° latitude-longi- 
tude grid. Multi-level data were interpolated to 17 stan- 
dard pressure levels by each modelling group, with the 
exception of parameterization forcing (PF) data, which 
were requested on model levels to retain the complete 
vertical structure. 

For the diagnostics presented here and in the ATLAS, 
data for individual models are plotted on the submitted 
grid without further interpolation. Multi-model means 
and standard deviations are obtained by first linearly 
interpolating to a common 1° latitude- longitude grid 
before averaging the available models. 

3. Participating Models 

Table 2 lists the 16 models that participated in the 


APE intercomparison. These include a number of 
established models used in production appl ications for 
both numerical weather prediction (NWP) and climate, 
including the Intergovernmental Panel on Cl imate 
Change (1PCC), and a number of novel models less 
tested and iterated in real-world applications. The APE 
intercomparison may therefore be expected to produce 
a wider spread in many circulation metrics than for 
previous intercomparisons such as AMIP. 

Table 2 includes the commonly accepted model 
name, the group that contributed the data and the 
location of the group's home institution: more details 
including names of the contributing modellers are 
provided on the APE website. The group symbol 
is used to identify the models in figures, tables and 
discussions . The final col umn identifies whether each 
model was optimised to give top of atmosphere radia- 
tive balance for present day Earth climate . This will be 
used in assessing the global energy budget for the APE 
CONTROL experiment. 

The partici pating models span a range of hori- 
zontal and vertical resolution reflecting their climate 
and NWP applications. In addition there is one global 
cloud-system resolving model, although this was only 
integrated for a 30 day period. The models include a 
diverse range of both dynamical core and parameter- 
ization choices. Table 3 provides some basic charac- 
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Table 3. Dynamical properties of palticipating models 


GROUP SYMBOL 

DYNAMICAL CORE 

WATER VAPOR 

HORIZONTAL 

VERTICAL 



TRANSPORT 

RESOLUTION 

RESOLUTION 

AGU 

Eulerian spectral 

Eulerian spectral 

T39 

L48 

CGAM 

lat-Jon grid point 

Eulerian grid 

3.75° x 2.5° 

L30 

CSIRO 

C-C 1 semi-Lag 2 

semi-Lag 

-210 km (C48) 

L18 

DWD 

icosahedral grid 

semi-Lag grid 

-JO 

L3 1 

ECM-CY29 

semi-Lag spectral 

semi-Lag grid 

T159 

L60 

ECM-CY32 

semi-Lag spectral 

semi-Lag grid 

T159 

L60 

FRCGC 

icosahedral Eulerian 

Eulerian 

-? km 

L54 

GFDL 

Jat-lon finite volume 

finite volume 

250 x 20 

L24 

GSFC 

lat-lon grid point 

Eulerian centered 

3.75° x 3° 

L34 

K1JAPAN 

Eulerian spectral 

semi-Lag grid 

T42 

L20 

LASG 

Eulerian spectral 

Eulerian grid 

R42 

L9 

MIT 

cubed sphere 

Eulerian grid 

-280 km 

L40 

MRI 

Eulerian spectral 

Eulerian spectral 

T42 

L30 

NCAR 

Eulerian spectral 

semi-Lag grid 

T42 

L26 

UKMO (N48) 

semi-Lag lat-lon grid 

semi-Lag 

3.75° x 2.5° 

L38 

UKMO (N96) 

semi-Lag lat-lon grid 

semi-Lag 

1.875° x 1.25° 

L38 


1 C-C denotes conformal cubic 
2 semi-Lag denotes semi-Lagrangian 


teristics of the models including the type of dynamical 
core, the method used for water vapor transport if it 
differs from that of the dry dynamics, the horizontal 
resolution and the number of model levels . Table 4 
lists major choices in the model parameterizations of 
convection and boundary layer turbulence. These can 
only provide labels for the different model character- 
istics that are indicative of their definitions and serve 
as a reminder. The actual schemes must be determined 
from the additional information and references on 
model formulation presented in the ATLAS . 

4. Forcing of the global circulation 

Before analysing circulation statistics from the aqua- 
planet simulations, the forcing of the aqua-planet circu- 
lation is compared with that of Earth. Although APE is 
configured as a perpetual equinox, it is compared with 
annual averages for Earth, to average over the observed 
large annual cycle of poleward energy transport and to 
avoid differences associated with seasonal uptake and 
release of heat by the ocean . Figure 1 shows the equi- 
noctial insolation for APE, which varies simply with 
cosine of latitude , and the annual average for Earth as 
calculated by the NCAR model. The APE insolation 
has a larger latitudinal gradient than the annual mean 
at all latitudes, so tends to force a stronger circulation. 
By prescribing SST, the fraction of insolation reaching 


the surface is lost, although an implied ocean pole - 
ward heat transport can be computed to balance the 
net surface energy flux. The SST profile for the APE 
CONTROL experiment, also shown in Fig. 1, is more 
peaked about the equator than the annual average for 
Earth, with unlimited mois ture availability at all longi- 
tudes: this is also expected to force a stronger meridi- 
onal circulation in the tropical atmosphere. In contrast, 
the APE SST is constant beyond 60° latitude, opposing 
the effect of increased insolation gradient on poleward 
energy transport there . 

Figure 2 shows the poleward energy transport for the 
APE CONTROL SST experiment, calculated following 
Gleckler et al. (1995). The multi-model mean is shown 
as solid lines and the standard deviation of the models 
about the mean is shaded. The total transport by the 
atmosphere and an implicit underlying ocean is that 
required to balance the net top of atmosphere (TOA) 
energy flux ; the oceanic transport is that required to 
balance the net surface energy flux and the atmospheric 
transport is the difference between the two. Globally 
averaged net fluxes are not relevant to the transport 
and have been removed: the global energy budget is 
discussed in Section 7. Observational estimates of 
the annual average poleward transport for Earth from 
Fasullo and Trenberth (2008) are included in Fig. 2 for 
companson. 
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Table 4. Parameterizations of participating models 


GROUP SYMBOL 

PBL 

SHALLOW CONVECTION 

DEEP CONVECTION 

AGU 

Mellor-Yamada 

None 

Emanuel 

CGAM 

Smith 

Gregory-Rowntrne 

Gregory-Rowntree 

CSIRO 

Holtslag-Boville 

None 

McGregor 

DWD 

Louis 

Tiedtke 

Tiedtke 

ECM-CY29 

Louis-Be ljaars 

Tiedtke 

Bechtold et al. 2004 

ECM-CY32 

Louis-Be ljaars 

Bechtold et al. 2008 

Bechtold et al. 2008 

FRCGC 

Mellor-Yamada 

None 

None 

GFDL 

Lock 

RAS 1 

RAS 

GSFC 

Louis 

RAS 

RAS 

KJ JAPAN 

Mellor-Yamada 

None 

Pan-Randall 

LASG 

Local vert diffusion 

None 

Manabe 

MIT 

Mellor-Yamada 

RAS 

RAS 

MRI 

Mellor-Yamada 

Randall-Pan 

Rmldall-Pan 

NCAR 

Holtslag-Boville 

Hack 

Zhang-Mcf arlane 

UKMO (N48) 

Lock/Richardson 

Gregory 1 990/Grant 

Gregory 1999 

UKMO (N96) 

Lock/Richardson 

Gregory 1990/Grant 

Gregory 1999 


1 RAS denotes relaxed Arakawa-Schubert. 
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Fig. 1. Incident solar short wave flux at the top 
of atmosphere (left axis), zonally averaged , 
comparing the equinoctia 1 insulation specified 
in APE (solid) and the annual average for Earth 
(dashed) . Sea surface temperature (right axis), 
zonally averaged , specified for the APE CON- 
TROL experiment (dash-dot) and the annual 
average for Emth (short-dash). Earth values m e 
averaged over northern and southern latitudes 
and SST is set to 1.8°C wider sea-ice. 


The total energy transport for the APE models is 
rather similar to Earth , reflecting similar TO A net flux 
distri butions, but the peak transport is slightly stronger 
and shifted poleward in APE. However, the differences 

uncertainties in the observational estimate (Fasullo and 
of approximately 0.5 PW are of the same order as 
Trenberth 2008) and, except at high latitudes, are less 
than the standard deviation of the APE models, which 
is large, reaching around 1 PW in mid-latitudes. The 
atmospheric transport dominates, with APE having a 

similar magnitude but an equatorward shift relative 
to Earth in mid-latitudes . In particular the APE trans- 
port is larger by around 1 PW in the tropics, giving a 
sub-tropical ’’shoulder” to the distribution. The atmo- 
spheric transport is significantly more constrained 
among the APE models than the total transport, with 
the inter-model standard deviation reaching only 0.4 
PW in the sub-tropics. This combines with the large 
spread in total energy transport to give a large standard 
deviation, around 1 PW, in the implied ocean transport 
required for energy balance. This implies that signifi- 
cantly different equilibrated climates would result if 
the APE models were coupled to either thermodynamic 
slab oceans or dynamic ocean models. 

The opposing latitudinal shifts in total and atmo- 
spheric transports are associated with significantly 
different net surface energy flux and implied ocean 
transport. For APE the ocean transport is weak in the 
tropics and increases through mid-latitudes, whereas 
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POLEWARD TRANSPORT 



Fig. 2. Poleward energy transport for the APE 
CONTROL experiment: multi-model mean (solid 
lines) and standard deviation (shading). Tota 1, 
ocean and atmosphere transports are inferred , 
respectively , from the top of atmosphere and 
surface net fluxes and their difference. Annual 
mean observational estimates for Ealth from 
Fasullo and Trenberth (2008) (dashed) . Values 
are zonal -time averages, averaged over northern 
and southern latitudes. 

for Earth it peaks in the sub-tropics and decreases 
through mid-latitudes. Despite the large spread, all the 
APE models share this distribution of ocean transport: 
most models have weak poleward transport throughout 
the tropics but three outlying models have significant 
equatorward transport and reduce the average to near 
zero (see the ATLAS). 

The differing latitudinal profiles of poleward trans- 
port for the APE CONTROL and Earth reflect different 
partitioning of the net TOA heat flux between heating 
of the atmosphere and ocean. At high latitudes the 
total poleward transport and its partition between 
atmosphere and ocean are insensitive to changes 
in the APE SST profile, as discussed by Williamson 
et al. (2013). The total poleward transport is -0.75 PW 
stronger in APE than for Earth at 60° latitude , forced 
by the stronger latitudinal gradient of insolation . The 
atmospheric transport is —0.75 PW weaker than Earth , 
presumably due to the absence of SST gradient pole- 
ward of 60° in APE. The warm isothermal SST also 
results in higher upward net heat flux at the polar 


surface, requiring larger oceanic heat flux convergence 
for energy balance and thus a much 1 arger poleward 
transport than Eaith, by -1 .5 PW at 60° latitude. 

In the deep tropics, the TOA net downward flux is 
around 75 W m- 2 for both APE and Earth. Of this, 
less than one third goes i nto the surface for the APE 
CONTROL, com pared with two thi rds for Earth (see 
Fasullo and Trenberth 2008, Fig.4 and the ATLAS). 
As the APE SST profile is flattened, going from the 
CONTROL to QOBS to FLAT cases defined by Neale 
and Hoskins (2000a), the TOA net flux changes little 
but the surface net flux increases substantial ly equa- 
torward of 15° and decreases poleward to around 40°. 
The tropical atmospheric transport therefore decreases 
and the implied ocean transport i ncreases with a rela- 
tively fixed total. The atmospheric dominance of trop- 
ical poleward energy transport in the APE CO TROL 
simulation is therefore ultimately determined by 
the more peaked SST profile than on Earth . This Is 
discussed further in Williamson et al. (2013). 

5. Zonal mean state 

In this section we consider the time average, zonal 
averages of the multi-model mean and the varia- 
tion among the models. The multi-model mean here 
includes all the models in Table 2 except FRCGC. 
Only a 30 day period was provided for the CO TROL 
experiment for this global cloud-system resolving 
model , due to its computational cost, so the model 
climate is unlikely to have reached equilibrium . The 
remaining models are equally weighted despite some 
models being outliers in various statistics . However, 
the mean includes a sufficient number of models that 
the impact of outliers is generally small. 

The variability among models is summarized by 
plots of inter-model standard deviation or is described 
in words without showing the individual models . 
These discussions are intended to emphasize the great 
variation among the models in most of the statistics 
examined . The corresponding plots are available in the 
ATLAS but specific figure numbers are not provided 
here . In general all the individual models differ signifi- 
cantly from the multi-model mean. 

The time average zonal averages are symmetric 
about the equator, which is to be expected given the 
hemispheric symmetry of the forcing. The two hemi- 
spheres are therefore averaged together to reduce the 
noise slightly and only one hemisphere is plotted , in 
the sense of the Northern Hemisphere . In addition , sine 
of latitude is used for the abscissa, to provide more 
detail of the tropical regions which exhi bit more struc- 
ture for many statistics. 
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5.1 Zonal wind and temperature 

The multi-model mean zonal average state shown i n 
Fig. 3 is very earth-like in structure, comparable to a 
symmetrized observational esti mate, but with signifi- 
cantly stronger westerly jets. The jet core is 62 m s- 1 
and occurs at 30° latitude and 200 hPa. The maximum 
low-level westerl y wi nd occurs only 5° poleward of the 
jet at 35°, so there is no separation between an upper 
tropospheric sub-tropical jet and the deep eddy-driven 
jet. This is due to the SST distri bution, which is more 
peaked about the equator than on Earth, as discussed i n 
Wi lliamson et al. (201 3), where the impact of changes 
in latitudinal SST profile on jet location and structure 
is considered. The zonal ly symmetric nature of the SST 
also contributes to the overall strength of the zonally 
averaged flow: the trade winds and low level polar 
easterly flow are also stronger than observed in the 
annual mean and in any individual season in ECMWF 
reanalyses (Kallberg et al. 2005). The polar easterl ies 
extend further equatorward towards the more equator- 
ward mid-latitude westerl ies. 

The inter-model standard deviation generally 
increases with height in the troposphere, along with 
mean wind speed. However, it is a relative minimum 
close to the jet axis, with peak val ues both poleward and 
equatorward of the jet co-located with the maximum 
wind gradients. This is generally due to a subtle pole- 
ward or equatorward shift of the jet in the individual 
models (seen as a dipole structure in differences from 
the multi-model mean in the ATLAS) rather than to 
a widening or narrowing of the jet. The peak in stan- 
dard-deviation is aligned vertically on the poleward 
flank of the jet but tilts poleward with height on the 
equatorward side, reminiscent of the stmcture of the 
leading mode of annular variability found in dynamical 
core simulations, which describes latitudinal displace- 
ments of the jet (see e.g., Sparrow et al. 2009). The 
standard deviation at the jet maximum is around 3 
m s- 1 • Individual model maximum values range from 
60 m s- 1 to 67 m s- 1 , with the exception of MIT which 
is the lowest at 54 m s- 1 • The peak standard deviation 
on the jet flanks reaches around 5 ms- 1 , while near the 
surface it is I m s- 1 or less. The standard deviation is 
somewhat higher on the equatorward flank of the jet , 
but the extent to which model variations in the Hadley 
circulation may contribute to this is unclear. 

In contrast to these modest values of inter-model 
standard deviation in the troposphere, the spread of 
zonal mean zonal wind increases significantly above 
150 hPa in the lower stratosphere. There is a local ized 
peak above the equatorial tropopause , but the largest 
values occur above the poleward flank of the jet. This 


is due to upward and poleward extension of the tropo- 
spheric jet into the stratosphere in several models, 
form ing a second zonal wi nd maximum near 50° at the 
model top, while a 1 arger number of models close off 
the tropospheric jet and form a well separated second 
maximum fort her poleward (evident in i ndividual 
model zonal wi nds in the ATLAS). This large varia- 
tion in stratospheric flow may be associated with the 
top boundary conditions of the model s. For example, 
the NCAR CAM3 applies a \1 2 horizontal diffusion 
at the top three levels of the model as a simple sponge 
to absorb vertically propagati ng planetary wave energy 
and also to control the strength of the stratospheric jets. 
The diffusion coefficient is chosen arbitrarily to yield 
reasonable polar night jets in earth-like si mulations 
(Collins et al. 2004). This diffusion provides a strong 
control on the jet and thermal stmcture near the top of 
the model. 

The multi-model mean zonal average temperature , 
shown in Fig. 3, is also superficial ly very earth-like. 
The temperature difference between the tropics and 
60° latitude is similar to the annual average on Earth, 
but the tropospheric temperature gradient is more 
limited in latitudinal extent and located further equa- 
torward , associated with the latitudinal SST profile . 
This is consistent with the stronger and more equa- 
torward westerly jet in the aqua-planet than in reality. 
Beyond 60° latitude the temperature gradient is weak 
in the aqua-planet by design, particularly in the lower 
troposphere . Here the polar atmosphere is less stati- 
cally stable than reality, so temperature gradients are 
closer to observed in the upper troposphere . 

The inter-model standard deviation of temperature 
generally increases with height, with a strong increase 
immediately below the tropopause at all latitudes. At 
each level there is a localized maximum in standard 
deviation , co-located with the maximum temperature 
gradient. In the upper troposphere this is close to the 
centre of the jet at 30° latitude and is co-located with the 
minimum standard deviation of zonal wind between the 
two peaks either side of the jet, consistent with thermal 
wind balance . The localized maximum extends down- 
ward following the maximum temperature gradient 
on the equatorward flank of the jet. Throughout the 
mid-troposphere the largest standard deviation is there- 
fore in the sub-tropics. This may also reflect the relative 
lack of thermal contact between descending air above 
the trade wind inversion and the underlying SST: it is 
notable that the standard deviation increases strongly 
between 925 hPa and 850 hPa in the sub-tropics. In 
contrast, equatorial deep convection communicates the 
SST throughout the troposphere . In mid-latitudes deep 
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convection also occurs in equatorward moving air in 
baroclinic waves, communicating the SST higher into 
the troposphere . Larger values of standard deviation 
occur over the polar cap, due to a variety of convective 
behaviour among the models there . 

The inter-model standard deviation of temperature is 
significantly larger in the stratosphere than the tropo- 
sphere, with a maximum value of 6 K in the region of 
strong veltical temperature gradient above the tropical 
tropopause. This is influenced by one outlier model 


but a number of models differ from the multi-model 
mean by almost 1 0 K (individual models are shown 
in the ATLAS) . Historically AGCMs have exhibited 
large cold biases at the polar tropopause, as discussed 
by e.g. Boer et al. (1 992) and Gates et al. (1999). There 
is evidence that use of a semi-Lagrangian dynamical 
core reduces the bias, although it remai ns a stubborn 
problem in some models (Williamson and Olson 1994; 
Chen and Bates 1996; Williamson and Olson 1998). 
Of the APE models, two of the three that use semi-La- 
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grangian dynamica 1 cores (ECMWF and UK.MO) are 
significantly warmer than the others i n the sub-polar 
lower stratosphere, on the poleward side of the tropo- 
spheric jet, confining the jet on its poleward upper 
flank. The CGAM model {HadAM3 with 30 levels) 
is also warmer than the multi-model mean in the 
pol ar lower stratosphere, consistent wi th the im pact 
of increased vertical resoluti on i n that model seen by 
Pope et al. (200 I ). The remaining semi-Lagrangian 
model , CSI RO, exhibits the opposite sense of tempera- 
ture anomalies in the lower stratosphere, possi bly due 
either to its lower vertical resol ution or i ts treatment of 
the vertical trajectory calculation. 

The majority of the features of the APE inter-model 
variation discussed above are also seen in AMJP si mu- 
lations of Earth climate performed as part of phase 3 of 
the Coupled Model Intercomparison Project (CMIP3, 
Meehl et al. 2007). These are plotted for the annual 
mean climatology in the ATLAS. The standard-de- 
viation of zonal wind scales approximately with the 
multi-model mean wind, which is lower in the AMIP 
annual mean than in APE by a factor of two. The zonal 
wind variation in AMIP is proportionally less than this 
only on the equatorward flank of the mid-latitudes jets, 
while it is higher in the tropical upper troposphere. 
Temperature variation among the AMIP models is 
generally similar to APE , but the peak variation above 
the tropical tropopause is somewhat smaller in AMIP, 
presumably reflecting the use of established and opti- 
mized AGCMs, compared to a mix of established and 
recently developed models in APE. The AMIP tempera- 
ture variation is larger than APE at low levels in the 
northern (but not southern) mid-latitudes , reflecting 
the interactive land surface in AMIP compared to 
prescribed SST in APE. 

5.2 Meridional circulation 

The Eulerian mean meridional circulation in Fig. 
3 shows a three-cell circulation in each hemisphere 
similar to that for Earth, with a strong Hadley cell in 
the tropics, a thermally indirect circulation that strad- 
dles the jet and surface westerlies near 35° latitude, 
further equatorward than for Earth, and a weak ther- 
mally direct polar circulation . The Hadley circulation 
is much stronger than for Earth at equinox, with merid- 
ional wind reaching 5 m s- 1 in the boundary layer and 
at 200 hPa for the multi-model mean . This is mainly 
due to the SST profile, as discussed earlier in Section 
4, but zonal symmetry of the SST and the proximity 
of the extra- tropical jet may also contribute. Strong 
upward motion is centered on the equator and down- 
ward motion extends from around 1 0° to just beyond 


30° latitude. The multi-model mean meridional wind 
shows a si ight hint of secondary cells in the vertical in 
the tropics between the main equatorward surface flow 
and the poleward upper tropospheric flow. There is 
also a barely discemable m id-tropopspheric minimum 
in the equatorial ascent. In fact all models except M IT 
show equatorward flow to some degree below the 
upper level poleward flow and many models also show 
a poleward flow above the low level equatorward flow. 

The standard deviations of mean meridional wind 
and vertical motion are superficially similar to their 
respective mean distributions, so variation among the 
models is primaril y due to differences in the strength 
of the meridional overturning . However there are also 
significant differences in both vertical and latitudinal 
structure. 

The deep structure of the standard deviation of mean 
meridional wind compared to i ts mean reflects varia - 
tions i n the vertical structure of the Hadley circulation 
among the models. The tropical upper tropospheric 
poleward flow peaks at 200 hPa in all but two models 
and the top of the equatorial upward motion varies from 
200 hPa to 150 hPa. However, the strength and depth 
of both the upper level outflow and the boundary layer 
inflow vary significantly. Variations in the strength of 
the secondary cells also contribute to the standard devi- 
ation of meridional wind near 700 hPa and 400 hPa and 
to the peak in variation of upper tropospheric vertical 
motion near 15° latitude . Here, some models show a 
distinct mid-level minimum descent, with a second 
peak in the upper troposphere. 

In the lower troposphere, the variation in low-level 
equatorward flow maximises equatorward of the multi- 
model mean, pointing to significant differences in the 
latitudinal structure of low level convergence in the 
JTCZ among the models . The variation in vertical 
velocity is largest in the equatorial region and is due 
to large variations in both the strength of the upward 
motion and its latitudinal structure. All models have 
ascent at the equator and the majority form a single 
maximum there, but the width of the ascent varies and 
several models form a double maximum straddling the 
equator. This differing structure is responsible for the 
very large variation at the equator itself and for the 
double peak in latitude, with minimum variation near 
50, 

In mid-latitudes, the maximum variations in merid- 
ional wind and vertical motion are shifted in latitude 
relative to their means. This is due to joint variations in 
the strength and latitude of the Ferrel cell: in a subset 
of models the jet is fmlher poleward and the Ferrel cell 
is stronger and extends fi.ul her poleward. 
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The different sub-grid scale parameterizations 
employed in the models are 1 ikel y to play a major rol e 
in dlivi ng the large variations i n the strength and spatial 
stmcture of the tropical meridional and vertical motion. 
However, sensitivi ty studies in the NCAR model by 
Wi 1 liamson and Olson (2003) and Wi 11 iamson (2008b) 
show that details of the i mplementation of the parame- 
terizations and the resolution at which they are applied 
may also lead to significant differences. In contrast , 
dynam ical processes play a more significant role i n 
variation of the mid-latitude Eulerian mean circulation 
(see Section 6.1 below). 

5. 3 Water vapour 

The multi-model mean specific humidity in Fig. 

3 is determined predominantly by the distribu- 
tion of temperature, but the narrow equatorial peak 
in mid-tropospheric humidity reflects the strong sub- 
tropical descent centered on 20° latitude in the APE 
CONTROL. Here the relative humidity is below 20% 
for the multi-model mean. On the equator, high 
relative humidity extends well above the boundary 
layer to 600 hPa. There is a separate maximum at the 
Hadley cell outflow level and a minimum below this at 
400 hPa . It is unclear whether this minimum is asso- 
ciated with the secondary cells discussed above, since 
the secondary inflow (equatorward meridional wind) 
is located below the minimum relative humidity. The 
height of the boundary layer is most strongly confined 
in the sub-tropics, by descent of dry air above . Further 
poleward , high relative humidity extends higher, due 
to boundary layer ventilation by baroclinic waves in 
mid-latitudes and by convection over the polar cap. 
Very dry air in the lower stratosphere, with relative 
humidity below 5% in the multi-model mean, extends 
above the tropopause temperature minimum at most 
latitudes . The exception to this is the polar cap, where 
colder air extends up to 10 hPa and the mean relative 
humidity exceeds 10%. 

In the troposphere the largest humidity variation 
among the models occurs at the sub-tropical boundary 
layer inversion and is seen in both absolute and relative 
humidity variation at 850 hPa . The boundary layer top 
is not well resolved in the interpolated isobaric data 
but , in a number of models 850 hPa is clearly above 
the boundary layer in relatively dry air, while in other 
models the relative humidity gradient occurs mainly 
between 850 hPa and 700 hPa (individual models are 
shown in the ATLAS). An additional level at 775 hPa 
would have enabled better resolution of these vertical 
stmctures . Relative humidity also varies widely among 
the models in the region of strong vertical gradient at 


the tropopau se. In the tropi cs the maximum variation 
is located at the cold point, I 00 hPa . This is distinct 
from the region of maxi mum temperature variation 
higher up, and is instead associated with valiations in 
the vertical extent of high tropospheric humidity. In 
the majority of models the relative humidity maximum 
is coincident with the poleward meridional flow 
maximum at or below 200 hPa, but i n three models 
(ECM-CY32, GFDL, K I JAPAN) it is located at 1 00 
hPa, above the layer of pol eward flow. This structure 
is worthy of farther i nvestigation. It coul d be associ- 
ated either with differences in the vertical profiles of 
convective heating and detrainment of moisture or 
with large scale advection in these models. There is a 
secondary maximum of humidity variation at 300 hPa, 
also associated with these surprisingl y large variations 
in vertical humidity structure in the equatorial upper 
troposphere. 

Globally, the largest variation i n relative hum idity 
among the models occurs in the polar lower strato- 
sphere at 70 hPa, where the mean relative humidity 
is only 11% but its standard deviation exceeds 40%. 
These statistics are significantly biased by three outlier 
models. MIT has large negative humidity there 1 and 
omitting it increases the multi-model mean to 21% and 
reduces the standard deviation to 24%. Two further 
models with greater than 50% relative humidity at 70 
hPa contribute the majority of this residual standard 
deviation . 

The regions of large humidity variation described 
above are associated with major differences in tropo- 
pause relative humidity structure in individual models 
(shown in the ATLAS) . At one extreme , the LASG 
hygropause is virtually flat between 150-200 hPa at all 
latitudes, due to this model 's very low vertical reso- 
1 ution (9 levels). In contrast , the ECMWF and MIT 
hygropause is above 100 hPa in the tropics but below 
200 hPa in the extratropics. At the poles, AGU is rela- 
tively dry in the upper troposphere but much moister 
above this at 1 00 hPa . 

The distribution of humidity variation among 
the APE models is generally similar to that for the 
AMIP annual average (shown in the ATLAS), taking 
account of the stronger Hadley circulation in APE and 
weaker tropical features in AMIP due to averaging 
over the annual cycle . Only two regions show signif- 
icant differences. First, in the equatorial upper tropo- 


1 LASG is the only other model with negative specific 
humidity (at the tropical tropopause) , although its relative 
humidity remains positive everywhere, suggesting interpola- 
tion problems from the 9 model levels. 
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sphere the region of large standard deviation extends 
further below the Hadley cell outflow in APE, despite 
the multi-model mean vertical profile being similar to 
AMIP. Second, the AMIP variation is much smaller in 
the polar lower stratosphere, presumably because the 
AMIP ensemble does not include such large outlier 
models , although circulation differences and the annual 
averaging may also contribute. 

5.4 Precipitation 

The single versus double equatorial ITCZ structure 
in vertical velocity is reflected in the zonal average 
total preci pitation shown in Fig. 4 for the individual 
models . The preferential forma t ion of a single or 
double equatorial preci pitation maximum has long 
been a feature of AGCM simulations of Earth climate 
and idealised aqua-planets, and it remains something 
of a modern modelling mystery. Due to feedbacks 
between low level winds and equatorial SSTs, an unre- 
alistic double-ITCZ is a major systematic bias in many 
coupled ocean-atmosphere simulations (see e.g., Lin 
(2007) and the discussion in Blackburn and Hoskins 
(2013)). In Fig. 4 the total precipitation is plotted equa- 
torward of 40° latitude and with the models split over 
two panels to better expose the tropical behaviours. The 
mid-latitude maximum occun-ing at around 35° latitude 
is associated with baroclinic waves . The curves in the 


unplotted region poleward of 40° show little variation 
between them as can be verified in the ATLAS. It has 
been verified that averaging the two hemi spheres does 
not misrepresent the ITCZ structures. All models wi th 
a double-ITCZ in the hemispheric average do contain a 
true double-ITCZ in the 3 -year mean . Using daily data 
for 1 year, the ITCZ structures appear stable down to 
monthl y timescales but relative north-south intensities 
do vary significantly on this timescale . 

There is a large varia tion in the maximum ITCZ 
preci pitation, from 10 to 34 mm d 1 • The CGAM, 
GSFC, NCAR and UKMO(N48,N96) models have 
clear double structures, where double is defined as the 
presence of a minimum equatorward of the maximum . 
The DWD , FRCGC and ECM-CY29 models have 
a single maximum at the grid point on the equator. 
This might be influenced by int erpolat ion, since these 
three models submitted data interpolated to uniform 
latitude-longitude grids. For the ECMWF model, an 
equatorial maximum represen ts a maximum at the first 
Gaussian latitude on each side of the equator. The A GU , 
GFDL, LASO, K1JAPAN, MIT and MRI models have 
no equatorial grid point , but form a single structure 
across the two points spanning the equator: thus the 
curves are flat approaching the equator and represent a 
single structure . The CSIRO and ECM-CY32 models 
form a broader single structure across the equator and 
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two adjacent points (one in each hemisphere), i.e. the 
equatorial value is not significantly different from the 
adjacent points . Again , this might be influenced by 
interpolation. 

Figure 5 shows the tropical-average total precipi- 
tation for the individual models, obtained by zonally 
and temporally averagi ng between 20 °S and 20°N, 
and the contributions from convective and large-scale 
processes. The tropical-average total precipitation is 
around 4.5 mm d 1 , which corresponds approximately 
to a radiative convective equi 1 ibrium averaged over the 
entire tropics 2 . The range among the models , 1 .2 mm 
d- 1 or 25-30% of the mean, is modest but significantly 
larger than the range for the global average shown 
below in Section 7. This is due to differences among 
the models in several processes: radiative cooling 

rates, the latitudinal extent of the Hadley circulation 
descent region that is in exact mass balance with the 
ITCZ ascent and transient eddy poleward moisture flux 
in the sub-tropics. Convective preci pitation dominates 
the total in all models, but the variation in large-scale 
precipitation fraction is large , ranging from negligi ble 
in the UK Met Office family of models to almost 50% 
in K 1JAPAN. FRCGC, with a horizontal resol ution of 
7 km, does not include a convective parameterization. 

The modest variation in tropical-average total precip- 
itation also contrasts with the large variation in preci pi- 
tation intensity at the equator, indicating compensating 
effects of ITCZ intensity and width. This is evident in 
the latitudinal profiles in Fig. 4, where precipitation 
curves for individual models cross close to 4 ° latitude, 
giving minimum spread there . Therefore, the large vari- 
ation in ITCZ structure and intensity among the models 
occurs despite a relatively strong budget constraint on 
tropical-average precipitation . 

6. Maintenance of zonal mean state 

This section considers the processes that maintain 
the zonal mean state, focusing on the transient eddy 
dynamical fluxes and components of the parameteriza- 
tion forcing. The primary interest is not the detailed 
budgets , which must balance closely for the 3 year 
climate averages in each model , but rather the spatial 
structures of the dominant processes and their variation 
among the models. The relevant diagnostics are avai 1- 
able for the majority of models, apart from the param- 
eterization forcing from individual processes . Details 
are given below. 


2 Latent heating from 4 nun d 1 balances a radiative flux 
divergence from the atmospheric column of around 120 W 
m 2 equivalent to a cooling of 1 K d 1 throughout the depth. 
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6. 1 fransient eddies 

The zonally avemged covariance of the quantities a 
and fJ is partitioned into stationary zonal mean (SM), 
transient zonal mean (TM), stationary eddy (SE) and 
transient eddy (TE) components, respectively, in 

[a.BJ = [a][,B]+[a]'[,B] r +^ *.£*y +[a' *B'*] 

where the overbar ( ) denotes the time average and 
prime (S the deviation from it, square brackets [ ] 
denote the zonal average and star ( *) the deviation 
from it. 

SM quantifies the zonal-time average state and the 
variances are simply the squares of variables in Fig. 3. 
TM quantifies time variations of the zonal average state 
and is small for many of the variables discussed below. 
For Earth both stationary eddy (SE) and transient eddy 
(TE) fluxes contribute to the mean state budgets . For the 
APE CONTROL, forced by zonally symmetric SST, 
SE contributions are small, so the zonally symmetric 
overturning and transient eddies dominate the trans- 
ports and budgets . The Hadley circulation has already 
been shown to be stronger than for Earth, consistent 
with the atmosphere contributing a greater fraction of 
total poleward energy transport in the tropics for APE 
than for Eaith (Fig. 2). In the extratropics the APE and 
Eailh atmospheric transports are more comparable , 
so the APE TE fluxes are expected to be of similar 
magnitude to the total eddy (SE + TE) fluxes for Earth. 
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The TE fluxes have not been time filtered, so synoptic 
ti mescale eddies and lower frequencies contribute to 
the diagnostic s in this section. Wavenumber-frequency 
anal ysi s of 250 hPa meridional wind at 30° latitude, 
shown in the ATLAS, confi rms that a wide range of 
ti mescales is present in the storm-tracks in al 1 the 
models. The 3 year samples do contain a zonal wave- 
number 5 structure that has non-zero SE transport, but 
the noise of a single 3 year sample is too large to allow 
systematic comparison of the models. This feature is 
discussed separately in Section 9. 

Figure 6 shows the multi-model mean and the inter- 
model standard deviation, to summarize the variabi lity 
among the models, for selected TE cova- 
riance quantities. The variances [C ll’ *) 2 j , [C V f *) 2 ], 

^ \ ■ J 
[( w' ) 2 J jshow the location and magnitude of transient 

eddy activity. The covariances [w'*v'*], [v'*T r *] 
and [V * q ' *] show the TE poleward fluxes of westerly 
momentum, temperature and moisture. The SM, TM 
and TE components of these and additional quantities, 
including vertical fluxes, are shown in the ATLAS. 
Table 5 1 ists the mid-latitude maximum value of each 
quantity for the eleven (out of sixteen) models that 
submitted covariance data. 

The multi-model mean TE covariances in Fig. 
6 are mostly dominated by transient activity in the mid- 
latitude storm-tracks and resemble those for Earth 
(Kallberg et al. 2005) . The kinetic energy (KE) compo- 
nents [( u 1 * ) 2 ] and [( v' * ) 2 1 peak in the upper-tro- 
posphere while [C w ’*) 2 ] peaks in mid- troposphere 

directly below [( u ) 2 ] . The magnitudes of all three 
quantities vary greatly among the models, with stan- 
dard deviations generally 20-25% of the multi-model 
mean . The maxi mum val ues, shown in Table 5, vary by 
roughly a factor of two, with standard deviations also 
in the range 20-25%. The GSFC values are generally 
the smallest, i ndicating least transient activity, whil e 
K I JAPAN is among the largest , wi th particularly large 
transient KE. Values do not vary systematically with 
horizontal resolution, implying that factors such as 
numerical damping and parameterization character- 
istics help to dete rmine eddy ampl itudes. The ratio 
[(v'*) 2 1/[( u' *) 2 j in Table 5 is greater than unity 
in all the models except DWD , but varies widel y 
between I . I 0 and 1 .52. This ratio appears (indirectly) 
in the E-vector, E = ( V 12 -U 12 , -u'v' ) , defined by 
Hoskins et al. ( 1983), so its large range implies that 
the horizontal anisotropy and zonal group propagation 
characteristics of transient eddies vary widel y among 
the models. This could be associated either with high 
frequency transients alone or with differences in the 


residual of opposing high and low frequency contri- 
butions to the zonal E-vector component, which is a 
feature of observations found by Hoskins et al. ( 1983). 

The poleward flux patterns for the multi-model mean 
in Fig. 6 are again sim ilar to those for Earth. Now the 
corresponding standard deviations (SDs) are signifi- 
cantly different from the means. For \u' *v' *] the mean 
and SD are both largest in the upper troposphere, but 
maximum SD coincides with the maximum gradient 
between the poleward and equatorward flux maxima. 
For the SD maximises on both sides of the 

upper tropospheric peak in the mean. This quadrature 
relationship between mean and SD reflects variations 
in both the magnitude and location of these particular 
TE fl ux structures among the models. For momentum 
flux, the standard deviation is again around 20% of the 
multi-model mean in Fig. 6, and also in Table 5 for 
the poleward and equatorward fl ux maxi ma. The pole - 
ward flux maximum at 30° latitude is stronger than the 
equatorward flux maximum at 50° in all models, but 
the ratio varies widely, between 1.2 and 2.2 in all but 
the MIT model , for which it is 4.0. This indicates large 
differences in meridional propagation of TE wave 
activity in the storm-tracks , in addition to the zonal 
variations discussed above. 

The variation among the models is considerably 
smaller for the lower tropospheric temperature and 
moisture fluxes, with the SD of the maximu m val ues 
in Table 5 being only 10% of the mean for [v 1 * T' *] 
and 13% for [v' *q' *] . This points to stronger budget 
constraints on these quantities than for momentum . 
There is scope for further analysis of the fluxes and 
vertically integrated budget constraints . 

A number of the TE covariances additionally show 
tropical structure . This is generally of smaller ampli- 
tude than in mid-latitudes , with the notable exception 
of vertical velocity variance, for which the multi- 
model mean in Fig. 6 is similar in mid-latitudes and at 
the equator. Moreover, the standard deviation is signifi- 
cantly larger than the mean at the equator, whereas 
it is only 20-25% of the mean in mid-latitudes . The 
large differences in mean ITCZ intensity and structure 
among the models, discussed previously, are therefore 
accompanied by large differences in the transient zonal 
asymmetries . 

The horizontal wind variances reveal that the differ- 
ences in vertical velocity variance are associated with 
winds in the equatorial upper troposphere, at the same 
level as the mean Hadley cel 1 outflow, 200 hPa. The 
variabil i! Y is domi nated by east-west winds, since 
r(w'*) 2 ] exceeds [(v'*)2] by a factor of around 
four (note the different colour scales in Fig. 6). This 
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suggests significant differences in the equatorial wave 
spectra, principally in low zonal wavenumbers and/or 
wave types for which zonal flow domi nates. Although 
LASO has significantly higher tropical variance than 
other models, it does not dominate the inter-model 
standard deviations. 

o — V-, — — 

The covariances [u’*v'*] and [v’*T'*] also show 
significant variation about the multi-model mean i n 
the tropical upper troposphere. The momentum flux 
is downgradient, with equatorward flux of westerly 


momentum towards the peak easterly wind at the 
equator. The mean temperature flux is small but also 
equatorward. In addition there is a poleward mois- 
ture flux maximum near 10° latitude throughout the 
mid-troposphere in the multi-model mean, with weak 
equatorward flux near the surface. The variation 
among the models is more confined to the lower tropo- 
sphere and peaks near 5° latitude, presumably associ- 
ated with the differences in zonal-time average ITCZ 
structure. These fluxes may be evidence of equatorial 
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Table 5. Maximum Transient Eddy statistics 


model 

te uu 

te w 

te omom 

te vt 1 


te uv 

tevq 

units 

m2 s-2 

m2 s-2 

x IO 3 Pa^ 

K m s- 1 


m2 s-2 

111 s 1 g kg I 

AGU 

330 

451 

48 

19 

78 

-39 

8.0 

CSIRO 

309 

438 

41 

20 

59 

-33 

8.8 

DWD 

416 

357 

54 

19 

50 

-31 

7.3 

ECM-CY29 

404 

444 

56 


65 

-42 

7.6 

ECM-CY32 

417 

459 

57 

19 

65 

-44 

7.3 

GFDL 

274 

330 

36 

16 

57 

-41 

6.7 

GSFC 

281 

302 

24 

15 

37 

-31 

5.6 

K1JAPAN 

453 

690 

51 

19 

73 

-37 

8.1 

LASG 

289 

382 

48 

21 

76 

-34 

6.9 

MIT 

264 

370 

53 

21 

61 

-15 

8.8 

NCAR 

370 

544 

40 

20 

63 

-44 

8.5 


transient eddy was not available forECM-CY29. 


wave activity in the presence of mean state gradients, 
not part of the classical theory by Matsuno and Gill 
which requires a resting basic state . The equatorial 
wave spectra are discussed in more detail in section 
1 0. I. The transient mean covariances (TM, shown in 
the ATLAS) reveal significant temporal variation of the 
zonally averaged tropical moisture and Hadley circu- 
lation, and that these contribute to the total moisture 
fluxes. 

6.2 Parameterized thermo dynamic. forcing 

Parameterization forcing diagnostics were not 
submitted for all the models and, of those that were, 
not all incl uded the complete set ofrequested variables. 
For this reason and to avoid loss of vertical structures 
by interpolation to a standard grid, a multi-model mean 
has not been computed. Instead, selected variables are 
shown for individual models on the model vertical grid. 
Only heating terms that are well defined are presented 
here. Forcing by moist processes is partitioned differ- 
ently between models and so individual parameteriza- 
tions cannot be sensibly compared . More complete 
forcing diagnostics are included in the ATLAS. 

Figure 7 shows zonal and time averaged total param- 
eterized heating in the simulations . The largest inter- 
model variation is in the tropics . Heating variations are 
much smaller in the mid-latitude storm-track region, 
except for AG U with about 25% stronger heati ng and 
GSFC with double the heating rate . 

Variations in equatorial heating are associated with 
the convective and cloud parametelizations (shown i n 
the ATLAS) and for individual models are consistent 
with the variations in ITCZ vertical motion and precip- 


itation seen earlier. Most of the models have minimum 
heating near 600-750 hPa, perhaps associated with 
evaporation below the freezing level or with outflow 
from convection that reaches mid-troposphere (repre- 
senting cumulus congestus). Each model 's heating 
minimum is below its mid-level minimum in vertical 
motion (cf. Fig. 3) and closer to the level of secondary 
poleward mean flow in the low/mid troposphere, to the 
extent that these dynamical structures are resolved by 
the standard isobaric level data. 

In the sub-tropical subsidence region , large varia- 
tions in the rate of cooling at the boundary layer inver- 
sion are dominated by longwave radiative cooling, 
shown in Fig. 8, partly offset by qualitatively similar 
but weaker variations in shortwave heating (shown in 
the ATLAS). This large variation in radiative cooling is 
associated with tropical low level clouds, which have 
been identified as important for estimates of climate 
sensitivity. Medeiros and Stevens (20 1 1) argue that 
in aqua-planet simulations these are shallow cumulus 
clouds , not stratocumulus , and that they occur in similar 
dynamical environments to those in Earth simulations 
and in observations. Therefore, the shallow cumulus 
regime in the APE simulations can be compared with 
observations, just as well as the earth-1 ike configura- 
tions can be. 

Above the sub-tropical inversion, the parameteriza- 
tion cooling in Fig. 7 is in better agreement between 
the models . However, its latitudinal structure is non- 
trivially determined, since longwave cool ing (Fig. 

8) peaks fi.l hher equatorward while shortwave heating 
mini mises further poleward. The total parameterized 
cooling does not extend as far poleward as the mean 
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descent or, to a lesser extent, the minimum in relative 
humidity in Fig. 3. This implies that cooling due to 
divergence of the transient eddy poleward temperature 
flux (Fig. 6) makes an increasing contribution to ther- 
modynamic balance with latitude in the sub-tropics. 

Two models include heating structures that differ 
significantly from the other models. In the boundary 
layer, parameterization heating is stronger and shal- 
lower i n K UP AN than in the other models shown i n 
Fig. 7. Turbulence is the only contributing process i n 
K 1JAPAN , while convection deepens the boundary 
layer heating in the other models (individual process 
contl ibutions are shown in the ATLAS). GSFC exhi bits 
a large ampl itude grid-scale structure i n the lower 
tropospheric heating profile in Fig. 7. This is due to the 
convection and cloud parameterizations, partly offset 
by similar structures in longwave radiation and turbu- 
lent m ixing. Sub-grid turbulence is also active above 
the equatorial boundary layer in AGU, opposing grid- 


scale heating structures due to convection and cloud. 

All models except LASG show longwave heating 
at the tropical tropopause in Fig. 8. LASG has very 
coarse vertical resolution in this region and probably 
cannot resolve the processes involved . The heating 
is due to absorption of upwelling radiation from the 
warmer troposphere, so is stronger in the sub-tropics, 
where clear skies and low tropospheric humidity lead 
to a larger upwelli ng longwave flux than from colder 
deep convective cloud tops at the equator. The heating 
valies significantly among the other models, despite 
both ozone and carbon dioxide concentrations being 
prescribed in the APE protocol. Thuburn and Craig 
(2002) found that, although longwave heati ng at the 
tropopause col d point is domi nated by ozone absorp- 
tion, CO 2 domi nates the sensitivity of heati ng with 
respect to temperature in the tropical substratosphere, 
between the convection top and cold point Whichever 
gas dominates the tropopause heating variation in the 






3 6 


Journal of the Meteorological Society of Japan 


Vol. 91A 


AGU 


DWD 


ECM-CY29 


ECM-CY32 



80 60 40 20 

LATITUDE 


0 80 60 40 20 

LATITUDE 


0 80 60 40 20 

LATITUDE 


f u r Tr I lt T 

0 80 60 40 20 0 

LATITUDE 


GSFC 


K1 JAPAN 


LASG 


NCAA 



1000 III I I I I 1 
80 60 40 20 

LATITUDE 


0 80 60 40 20 

LATITUDE 


u i i i — i — i — r 

0 80 60 40 20 

LATITUDE 


80 60 40 20 

LATITUDE 


UKMO(N96} 



_*r v ,v/ r ?? 

80 60 40 20 

LATITUDE 


— 

0 80 60 40 20 

LATITUDE 
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APE models, Fig. 9 shows that, to first order, the long- 
wave heating variations damp differences in tempera- 
ture between the models . Shortwave variations (shown 
in the ATLAS) are small, so the longwave variations 
appear in the total parameterized heating (Fig. 7). This 


implies that the variations in tropopause temperature 
structure are dynamically maintained , by differing 
numerical computations of adiabatic cooling associ- 
ated with ascent. There is no obvious association of 
heating with veltical resolution in Fig. 9, but this inter- 
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pretation is consistent with the finding by Williamson 
et al. ( 1998) that tropical tropopause temperature and 
stmcture are highly sensitive to dynamical formula- 
tion at vertical resolutions typical of those in the APE 
models. 

7. Global average budgets 
7.1 1 Water 

Figure 1 0 compares the global average water budget 
displayed as bar plots for each model, calculated 
from the daily average, global average time series. 
The multi-model mean precipitation is almost 3 mm 
d- 1 , very close to observational estimates for Earth. 
There is a 20% variation in total preci pitation among 
the models, smaller than the 25-30% range of trop- 
ical averages seen in Fig. 5. A much larger variation is 
evident in the partition between convective and large 
scale stable preci pitation. Neglecting FRCGC , which 
has no convective parameterization , convective precip- 
itation varies by a factor greater than two and large scale 
varies by a factor of three. The temporal variability of 
global average precipitation is shown by the whiskers 
in Fig. 10, which range from plus one to minus one 
standard deviation of the time series of daily average 
values. The standard deviation is around 4-5% in most 
models, with the exception of DWD which shows very 
1 ittle day-to-day variation i n all global averages, most 
likely due to a loss of precision in the output data. 


Most models conserve water over the three year 
integration period, with evaporation minus precipi- 
tation being zero or nearly so. Exceptions are DWD, 
ECMWF and LASO, whose water budget residuals are 
approximately 1 .3%, 2% and 4% respectively of global 
preci pitation or evaporation . The main contributor to 
the ECMWF residuals is believed to be interpolation to 
a regular latitude-longitude grid before global integra- 
tion, since climate simulations of the ECMWF model 
produce a residual of about -0.4% on the model grid, 
much closer to those of the remaining APE models . 
A loss of numerical precision in the post-processing 
is believed to be the main contributor to the DWD 
residual. Advection and adjustment of negative vapour 
account for the LASO residual, the largest at -4% of 
global precipitation. The 30 day period for FRCOC is 
too short to rule out a change in the water vapor content 
of the atmosphere dominating the average. OFDL, 
OSFC, MIT and NCAR are the closest to zero. We note 
that NCAR has a fixer to ensure water is conserved . 
The remaining models show a minimal residual value . 

7.2 Energy 

Prescribing both insolation and SST in APE leads to 
an unknown global imbalance in TO A and surface net 
energy fluxes, even for models optimised for simula- 
tion of Earth cl imate. Atmospheric energy storage is 
negligible in multi-annual averages, so the imbalance 
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would heat or cool an underlying ocean. On Earth the 
present-day global imbalance is estimated to be 0.9 
W m 2 downwards , causing global warming (Hansen 
et al. 2005; Trenberth et al 2009), but it may be much 
larger for each of the aqua-planet SST distributions . 

The TOA and surface residual energy fluxes are 
shown for each model in Fig. II, together with the 
cloud fraction and albedo. For the multi-model mean , 
the TOA residual flux is around 13 W m 2 downwards, 
going into the ocean. The range across the models is 
surprisingly large, greater than 30 W m- 2 , even omitting 
FRCGC which is not in equili brium, and despite the 
simplified aqua-planet configuration which might be 
expected to constrain differences between the models. 
Several factors may contribute to the large range. 

First, experimental and numerical artifacts are rela- 
tively minor contributors. Most of the models accu- 
rately implemented the prescribed global insolation, to 
within 0.1 W m- 2 , but a small number did not, either 
by prescribing a different value of the solar constant 
or by using perpetual vernal equinox for Earth’s ellip- 
tical orbit 3 . The largest discrepancy is less than 4 W 
m 2 in global insolation, contributing only around 3 


3 We strongly recommend that future APE simulations use 
a circu lar orbit with zero obliquity. The Earth-Sun distance 
at vernal equinox is less than the annual average, increasing 
solar irradiance relative to the prescribed "solar constant". It 
is also difficult to identify the precise point in the calendar at 
which insulation is symmetric about the equator. 


W m 2 to net downward shortwave flux after al lowing 
for global albedo. Lack of formal energy conservation 
by the model atmospheres leads to a range of atmo- 
spheric energy gain/loss of 5 W m 2 , calculated as 
the difference between the TOA and surface residual 
fluxes . This contributes little to the range of residual 
flux, which differs by only 1.5 W m- 2 at the TOA and 
surface, although it contributes more to differences 
between individual models. 

Second, only a subset of the models listed in Table 
2 has been optimised to give global energy balance for 
present day Earth climate. The range of TOA residual 
flux is only moderately reduced for the optimised 
models but, omitting DWD which was optimised for 
short N WP timescales, the range is approximately 
halved to around 1 5 W m 2 . However, the contribution 
of this effect is uncertain , because the global imbal- 
ance for Earth climate in the non-optimised models is 
unknown . 

Third, optimisation of TOA flux for Earth climate 
may compensate differing systematic biases over land 
and ocean. Removi ng land in the aqua-planet then 
reveals the uncompensated bias over ocean, which will 
differ between models . T o quantify this effect, it would 
be necessary to compute fluxes averaged separately 
over land and ocean in AMIP si mulations for each 
of the models, and compare with the global residual 
fluxes in APE. 

Fourth, the APE CO TROL may be considered an 
extreme SST sensitivity experiment, with ealth-like 
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(rflux sfce, positive downwards) W m , cloud fraction (cld frac) and albedo. 


SSTs and land surface temperatures replaced by the 
specified aqua-planet values. This is analogous to the 
Cess type experiment (Cess et al. 1989, 1990). If the 
APE CONTROL SST is eliciting a range of sensitiv- 
ities in the global energy balance, we would expect 
the global imbalance to vary significantly for the 
gross variation of SST from the PEAKED through 
CONTROL and QOBS to FLAT aqua-planet prescrip- 
tions. This is certainly not the case for the multi-model 
mean (shown in the ATLAS), for which TOA residual 
flux varies by only around 2 W m- 2 . For individual 
models, over half approximately follow this mean vari- 


ation, though with varying magnitudes , while a smaller 
number show a reversed or less monotonic trend . The 
maximum variation with SST is 15 W m~ 2 (for AGU) 
but it is more typically around 5 W n> 2 , including the 
subset of models whose TOA balance has been explic- 
itly optimised for Earth cli mate. The inter-model range 
of residual flux is therefore rather insensitive to SST. 
It is si ightly larger for the PEAKED and FLAT cases, 
whose tropical climates are least si milar to Earth, as 
discussed in Williamson et al. (2013). 

The relative insensitivity of TOA global imbalance 
to aqua-planet SST suggests that the large range among 
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the models is mainly due to a lack of optimisation for 
Earth cl imate in the majority of the models and/or 
removal of differing biases in the TOA energy balance 
over land in the optim ised models. A combination of 
APE and AM IP si mulations by individual models 
would be required to quantify these effects and it is 
hoped that the inclusion of both AM IP and APE exper- 
iments in the protocol for the fifth phase of the Coupled 
Model Intercomparison Project (CMIP5, Taylor et al. 
2009) wi 11 facilitate this. 

Global cloud fraction and albedo 4 , also shown in Fig. 
I I , reveal that the variations in global energy balance 
among the models are dom inated by short wave reflec - 
tion by clouds. Variations in albedo correlate wel 1 with 
the TOA residual flux, and correspond to a range of 
46 W m- 2 in reflected (and net) shortwave flux. This 
is more than double the 20 W m- 2 range in outgoing 
longwave radiation (OLR, shown in the ATLAS), 
which is uncorrelated with the residual flux. However, 
variations in OLR tend to oppose those in the short- 
wave, leading to the smaller range in residual flux 
noted earlier, around 30 W m 2 compared to that in 
net shortwave. 

Global cloud fraction varies by a factor of almost 
two, from a particularly high value in AGU (0.87) to 
less than 0.55 in a cluster of four models (MRI and the 
3 UK Met Office models, CGAM and UKMO, again 
neglecting FRCGC). The latter group all have close to 
average albedo, so there is evidently a compensation 
between cloud fraction and cloud optical properties in 
this group relative to the remaining models . 

Figure 12 shows the latitudinal variation of zonal 
average cloud fraction and albedo from the individual 
models. Cloud fraction varies by more than a factor 
of two in both the sub-tropical and polar regions and 
is one of the most variable of all the metrics studied. 
In the sub-tropics, most of the models do cluster 
more closely, with 3 outliers having cloud fraction 
below 0.4 (FRCGC plus two from the group with low 
global cloud fraction), plus two in which the equato- 
rial maximum extends further into the sub-tropics. 
There is less general agreement at the poles. As was 
the case for the global averages, albedo is in some- 
what better agreement than cloud fraction, because of 
compensating cloud optical properties in a number of 
outlying models. The albedos in many models cluster 
within O.l of each other over most latitudes, but the 
range increases to more than 0.3 at the equator. This 
exceeds the fractional spread in equatorial OLR, so 


4 Global albe do is defined here as the ratio of the upward 
and downward global average shmtwave fluxes. 


shortwave variations dominate the energy balance at 
TOA and also at the surface there, drivi ng variations in 
pol eward energy tran sport i n the tropics , as discu ssed 
in Section 4. 

Cloud radiative forci ng (CRF) diagnostics were not 
collected for the APE models, so it is not possi ble to 
formally attri bute the differences in the TOA radiative 
balance among the A PE models to cloud versus clear 
sky effects. However, previou s studies have shown 
that the radiative response of clouds, and in particular 
sub-tropical boundary layer clouds, i s a major source 
of uncertainty in estimates of climate sensitivi ty in 
models (Bony and Dufresne 2005; Webb et al. 2006). 
Low -level clouds are expected to have a large radia - 
tive impact because their shortwave CRF dominates , 
whereas shortwave and longwave CRF more strongly 
oppose each other for deep tropical cloud. Variations 
in sub-tropical cloud fraction and optical properties 
among the APE models might therefore be expected to 
have the largest impact on the global energy balance. 
However this is not the case, since the TOA residual 
flux (shown in the Atlas) varies more in the ITCZ than 
in the sub-tropics. This large equatorial variation is due 
to the combination of the large range of albedo and the 
dominance of shortwave over longwave effects of deep 
cloud in the equatorial ITCZ noted above. 

8. Kinetic energy spectra 

Figure 13 shows the kinetic energy (KE) spectra of 
the models at 250 hPa , split over four panels for clarity. 
The figure includes the total kinetic energy (solid lines) 
and the divergent component (dashed lines), as a func- 
tion of 2-dimensional total spherical wavenumber n. 
The zonal average flow, i.e., zonal wavenumber m = 
0, is included in the KE for all n. Energy at a given 
total wavenumber n is associated with the spherical 
harmonics of vorticity and divergence at that wave- 
number . 

The KE spectra provide some indication of the 
damping from the dynamical core. The damping is 
due to explicit horizontal diffusion terms added to 
the dynamical equations, is inherent in the numerical 
approximations themselves , or is a combination of the 
two. In contrast to Earth simulations, which include 
large amplitude planetary waves forced by land-sea 
contrast and orography, the idealised APE configura- 
tion is dominated by a zonal average flow and transient 
eddies. It is hoped that this simpler circulation will help 
to isolate differences between models in the behaviour 
of the internal energy cascade process . 

The KE spectra in Fig. 13 are in good agreement 
among the model s for wavenumbers up to n - 15-20 
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Fig. 12. Zonal-time average cloud fraction (cld frac) and albedo. 


(apart from some exceptions discussed below) . At 
higher wavenumbers the effects of tranc ation and 
damping lead to a wide range of energies and spec- 
tral slopes. There is an approximate correspondence 
between the spectral slope in Fig. 13 and the transient 
eddy KE maxima in Table 5, discussed in Section 6. 1, 
with shal low spectra 1 slope corresponding to high 
transient eddy KE and v ice versa . Resolution i s an 
additional factor, with higher resolution contribu ting 
energy at additional wavenumbers, more so for spectra 


wi th shallow slope. Thus, the overall level of transient 
eddy KE is strongly controlled by the explicit and 
implicit damping in the dynamical core. 

The Eulerian spectra 1 transform models (AGU , 
K1JAPAN , LASO, MRI and NCAR) all maintain a 
slope of -3 to the truncation limit for the tot al kinetic 
energy. Generally, this is by design and obtained by 
adding a \7 4 diffusion term with the coefficient chosen 
to yield a slope of -3, to be consistent with observa- 
tional evidence of Nastrom and Gage (1985) at larger 
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length scales (n 80). In these cases with a -3 kinetic 
energy slope, the divergent component bas a slope 
of -5/3. It bas been argued that a \7 6 or \7 8 diffu- 
sion term which yields a steeper slope is desirable to 
reduce spectral ringing (MacVean, 1983) but \7 4 bas 
remained the most popular choice for Eulerian spec- 
tral transform models. The spectra from ECM-CY29, 
which is based on a semi-Lagrangian spectral trans- 
form core, drop off faster than -3 for the kinetic energy 
and faster than -5/3 for the divergent component. 
This behaviour is strongly influenced by the trunca- 
tion error of the (non-linear) semi-Lagrangian advec- 
tion operator and explicit or implicit diffusion in the 


dynamical core. The grid point models also all show 
faster drop-offs, some models faster than others, due to 
damping inherent in the grid point methods . The extra 
steep segment above wavenumber 40 of MIT might 
be due to the postprocessing interpolations to the lati- 
tude-longitude grid . Interpolation might also affect the 
spectra of CSIRO and DWD, but it bas littl e impact on 
the plotted component of the spectra for ECM-CY29 , 
i.e. for wavenumbers n < 90 that are resolved by the 
regular 2° interpolation grid. 

The model KE spectra are in much better agree- 
ment below wavenumber n 10. This is mainly due 
to similarity of the zonal average flow among the 
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models, which contributes most of the energy at low 
wavenumbers in the absence of significant stationary 
wave ampl itudes. Separate KE spectra for the zonal 
average state (zonal wavenumber m = 0) and eddy (m 
> 0) components in the NCAR model confirm this (not 
shown). 

The low wavenumber behaviour is consistent with 
a high resol ution aqua-pl anet simulation reported by 
Hamilton et al. (2008), hereafter HTO, (their Fig. 1 2). 
Compared to Earth simulations and observational anal- 
yses, also shown by HTO, the shape of the spectmm 
differs in APE, with lower energy in wavenumbers 2- 
4 and higher energy in wavenumbers I and 5-7 . Lower 
energy is consistent with weak stationary wave 
amplitudes, as concluded by HTO, but higher energy 
and the differing shape is associated with the latitu- 
dinal stmcture of the zonal average flow. In the APE 
CONTROL the westerly jet is stronger and at lower 
latitude than Earth, implying a stronger super-rotation 
and therefore higher energy in n — 1 (HTO, Fig. 12), 
and possibly higher energy at other wavenumbers. 
There is a large range of energy, around 5-30 m 2 s- 2 ,for 
a single wavenumber, n — 3, among the APE models 
in Fig. 13. Comparison with zonal-time average zonal 
wind for individual models in the ATLAS shows that 
a more poleward jet and weaker westerlies throughout 
the tropics are associated with higher KE in n = 3, and 
vice versa. Variations of the zonal average state among 
the APE models therefore project particularl y strongly 
onto the m = 0, n = 3 vorticity harmonic. 

Compared to Earth simulations, there is an exagger- 
ated difference between the energy in odd and even 
wavenumbers below n = 10 in APE, also seen in dry 
dynamical core simulations shown by Takahashi et al. 
(2006) and HTO. This is due to weak stationary waves 
in the idealised configurations and strong north-south 
symmetry of the zonal average state about the equator, 
which increases odd n and reduces even n vorticity 
harmonics for zonal wavenumber m = 0 (confirmed in 
the NCAR model, not shown). 

The divergent spectra in F ig. 13 show less agreement 
among the models at low wavenumber. Variations 
appear to be due to both the zonally symmetric over- 
turning circulation and low zonal wavenumber tropical 
variabi 1 ity (discussed later in Section 10). The Hadley 
cell outflow is generally above 250 hPa , so differences 
in both vertical structure and strength of meridional 
flow affect divergent KE at 250 hPa . LASO , which has 
the most active tropical variability, exhi bits the largest 
divergent KE at low wavenumbers. 


9. Mid-latitude low frequency variability 

A very noticeable wavenumber-five pattern often 
appears in mid-latitudes in long time averages of 
AP E simulations. It is most noticeabl e i n maps of the 
meridional velocity where wavenumber-five can be 
seen as the domi nant pattern. Examples are included 
in the ATLAS. Watanabe (2005), Cash et al. (2007) 
and Watanabe (2007) have discussed this mode and 
its relationship to the annul ar mode. Watanabe (2005) 
descri bes it as a slowl y propagati ng, wavenumber -five 
disturbance forced by high-frequency eddies. More 
recently, Zappa et al. (201 1) provide evidence that the 
feature is a weakly unstable barocl inic wave with very 
low phase speed, belongi ng to a barocl inic spectrum 
that obeys a well defined dispersion relation. I n the 
ECHAM5 model studied by Zappa et al (not incl uded 
in APE), the inverse energy cascade and phase locki ng 
with tropical convection each provide only a supple- 
mentary positive feedback. 

Unfortunately the design of the APE experimen t and 
the data collected do not allow a thorough examination 
of the phenomenon. As is seen in Fig. 14, the 3 -year 
averages of 250 hPa meridional velocity at 30° latitude 
show amplitudes of order 1 m s- 1 in zonal wavenum- 
ber- five. This gives an apparent amplitude of the mode 
at a single timescale, averaging over different phases 
due to slow phase propagation, the speed of which may 
differ between models . The true amplitude, fol lowing 
the phase propagation, cannot be robustly deduced 
from the data collected in APE, but is expected to 
decrease with increasing averaging period . The one 
year of 6-hourly data of a few two-dimensional fields 
is too short and does not allow examination of vertical 
aspects. However, the 3 -year averages do provide an 
estimate of the mode’s linear and covariance struc- 
ture in the APE models . The APE boundary condi- 
tions mean that there is nothing to set the phase of the 
wave in long-period averages and, indeed, apparently 
random phase variation is seen among the APE models 
and also for many 3 -year samples from the NCAR 
model (not shown). 

Nevertheless, most of the APE models indicate 
a significant mid-latitude wavenumber-five struc- 
ture in the 3 -year average which explains most of 
the ’’stationary” eddy variance on this timescale . The 
top row of Fig. 14 shows the amplitude of wavenum- 
ber- five of 3 -year average meridional velocity at 30°N 
and 30°S at 250 hPa. The left panel shows all the APE 
models except the very high resolution FRCGC which 
was only integrated for 30 days. The models are identi- 
fied numerical ly in Table 2. Clearly, there is no consis- 
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tency in the apparent amplitude of the mode, although 
this may be due to differences in both phase propagation 
and true amplitude among the models. The right panel 
shows the amplitude for 24 successive 3 -year averages 
from a 72 year integration of the NCAR model. This 
provides an indication of the variabil ity in ampl itude 
for 3 -year averages. The APE models are consistent 
with this va liability, with models 1, 8 and 12 slightly 
low in amplitude. The bottom row shows the percent 
variance explained by wavenumber- five at 30°N in the 


APE models (left) and in the 24 NCAR samples (right). 
The cases with relatively small amplitude show less 
variance explained. In the cases with more than 80% 
variance explained , maps of the meridional velocity 
shown in the ATLAS are clearly dominated by wave- 
number-five. In those around 60% (models 6, 14, 15), a 
strong wavenumber- five is modulated by other waves. 
In those below 40%, wavenumber- five is not visible. 

As found by Zappa et al. (2011) for the ECHAMS 
model , a number of the APE models also show weak 
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wavenumber- five modulation of the 3 -year average 
tropical circulation, including low level winds, evap- 
oration and moisture content in the sub-tropics and, 

less clearly, modulation of ITCZ precipitation (not 
shown). To see whether this produces a tendency for 

phase locking of wavenumber-five between northern 
and southern mid-latihldes, the top panels of Fig. 
14 show the amplitude of 3 -year average wavenum- 
ber- five at 30° and 30°S as a function of its harmonic 

phase difference between the hemispheres. For the 

CAR model the phase differences are clustered 
around 1 80° (meridional velocity out of phase, geopo- 
tential in-phase), suggesting significant phase interac- 
tion between the hemispheres. There is little evidence 
of phase interaction across all the APE models, but 

CAR does show higher apparent amplitude of wave- 
number-five than many of the other models. A lack 
of symmetry about 1 80° also suggests that the APE 
sample size is insufficient. 

Selected "stationary" eddy covariances for the 3- 
year averages are shown in the ATLAS for indi- 
vidual models, for all zonal wavenumbers and for 
wavenumber- five separately. For wavenumber- five the 
models generally have very weak poleward tempera- 
ture flux in the lower troposphere between -30-40° 
latitude, consistent with weak baroclinic instability . 
In the upper troposphere there is very weak poleward 
eddy momentum flux near 30-40° latitude, with equa- 
torward flux positioned further poleward in a few 
models. These structures are broadly similar to those 
of the dominant transient eddy fl uxes in Fig. 6. The 
wavenumber-five "stationary" eddy fluxes therefore 
appear consistent with the broader spectrum of baro- 
cl inic waves in the model storm-tracks, consistent with 
Zappa et al. (20 II), but anal ysis of longer integrations 
would be necessary for a more complete, quantitative 
analysis. 

IO. Tropical variability 

The overall magnitude of tropical variability in the 
models is compactly summarised by the standard devi- 
ation of six hour average preci pitation on the model/ 
data grid about the zonal-ti me average, shown in Fig. 
1 5. This varies by a factor of four among the models , 
excluding FRCGC whose standard deviation on its 
fi ne-scale grid is almost double the maxi mum among 
the remaining coarser scale models. The range of 
standard deviation , and the val ue i n many individual 
models , approximately scales with the zonal average 
precipitation at the equator in the ITCZ, which varies 
by a factor of around 3.5 in Fig. 4. Since there is a 
compensation between ITCZ intensity and width, due 


tpp n 



Fig. 1 5. Standard deviation of precipitation about 
the zonal-temporal average from 20° S to 20°N 
(tppn), mm day- 1 - The FRCGC value has been 
divided by 2 to show the variation among the 
models more clearly. 

to the tropical average energy budget discussed in 
Section 5.4, models with a narrow single ITCZ tend 
to have stronger preci pitation variabil ity, while models 
with a broader or double ITCZ tend to have weaker 
variabil i ty. This relationship bas been seen in a number 
of previous studies (e.g., Slingo et al. 1996; Wang and 
Schlesinger 1999; Lee et al. 2003, 2008). 

Differi ng horizontal resol ution among the models 
is likely to be a contributing factor to the i ntensity of 
preci pitation on the model grid, as Wi 11 iamson (2008a) 
found i n the NCAR model. The smallest horizontal 
scales oflatent heating are not directly relevant to large 
scale forcing, but they are part of the forcing spectrum 
whose interaction with the dynamics determines the 
nature of tbe tropical variabil ity, in particular its organ- 
isation in large scale equatorial waves . The following 
sections address these aspects of the simulations . 

10.1 Tropical waves 

The models show a particularly large vanatlOn in 
their equatorial wave propagation characteristics. This 
is most readily visualised in Hovmol ler plots of equa- 
torial preci pitation , averaged between 5°S and 5°N 
for sample 30 day periods for the individual models, 
in Fig. 16. It would be an understatement to say that 
the models display a wide variety of propagation char- 
acteristics! Eastward propagating features dominate in 
many of tbe models , with a phase speed that is remark- 
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ably constant among the models, around 15m s- 1 , and 
independent of zonal scale. Westward propagation 
dominates in only two models, AGU and MIT. lo a 
number of models, westward propagation of smaller 
scale features occurs within an eastward propagating 
envelope. Intense small scale quasi- stationary features 
occur in two further models, ECM-CY29 and GFDL. 
The intensity of the propagating features also varies 


widely. A diurnal cycle is visible in several models . 

This widely varying behaviour may be quantified 
using wavenumber-frequency diagrams, following the 
spectral analysis methodology of Wheeler and KiJadis 
( 1999). Figure 17 shows the log of the power of the 
symmetric component of the unnormalized spectra of 
the 6-hour averaged preci pitation averaged from 10° S 
to L0°N. The full power is plotted without removing a 






48 


Journal of the Meteorological Society of Japan 


Vol. 9 1A 


background spectrum in order to allow a comparison 
of the overall power of the waves. Nmmalization by 
a smooth background field, such as done by Wheeler 
and Kiladis, isolates spectral peaks that are often asso- 
ci ated with specific norma 1 modes or waves. However, 
si nee the normalization is done individually for each 
model, the signal associated with the overal 1 power is 
contained in the background spectrum and would be 
lost. In addition, i n the aqua-pl anet the wave character- 
istics emerge from the total without needing a normal- 
ization. Figure 1 7 i ncludes the conventional dispersion 
curves (labelled in the template panel) for odd merid- 
ional mode-numbered equatorial waves for equivalent 
depths of 12, 25, and 50m (see Wheeler and Kiladis 
1 999). The curves in the upper half of the plots show 
westward and eastward propagating inertia-gravity 
waves. Those at lower frequency show eastward 
propagating Kelvin modes and westward propagating 
equatorial Rossby modes. The period decreases with 
increasing equivalent depth in each set. The ordering of 
the models in Fig. 17 is the same as for the Hovmoller 
diagrams, for direct comparison. 

There are spectral peaks associated with the equa- 
torial wave modes in all the models, but the fraction 
of total power in these modes varies greatly. In a few 
models the modes are embedded in a background 
spectrum whose power increases with period, while 
in others most of the power projects onto the Kelvin 
and equatorial Rossby modes. A stronger background 
spectrum is characterised by more intermittent large 
scale eastward propagation in the Hovmoller diagrams, 
so occurs in those models which contain short lived, 
smaller scale westward propagation within an eastward 
propagating envelope. 

Eastward propagating Kelvin modes are present in 
all models , and dominate the spectrum in almost half. 
Kelvin mode power generally peaks at low frequen- 
cies but extends further to high frequencies in some 
models , while 2-5 day periods dominate in the NCAR 
model. These differences appear in the Hovmoller 
diagrams as variations in the zonal wavelength of east- 
ward propagating features whose phase speed varies 
remarkably little between the models. The small phase 
speed variations that do exist, e.g. faster in GSFC and 
NCAR, are visible in both figures. The agreement 
on eastward propagating phase speed is in marked 
contrast to aqua-planet simulations by Lee et al. (2008) 
in which AGCMs were coupled to a mixed layer ocean. 
Lee et al. 's finding of significant differences in phase 
speed may be due to widely varying equatorial SSTs 
among their coupled models, leading to differences i n 
static stability, precipitation rates and heating profiles. 


Westward propagation dominates in a smal lemumber 
of models, either as a westward bias in the background 
spectrum in LASG and MIT, or as a stronger projec- 
tion onto the equatorial Rossby modes in models such 
as AG U. However, westward low frequency power is 
generally less dispersive than the Rossby modes, with 
higher zonal wavenumbers '’lifted” to higher frequency. 
This is 1 ikely to be a doppler shifting due to the pres- 
ence of easterly zonal flow throughout the equatorial 
troposphere in the A PE CONTROL (Gui-Ying Yang, 
201 1, personal communication). 

A few models contain some clearly isolated waves, 
such as NCAR and to a lesser extent MRI with a I 0 
day eastward propagating wavenumber 6. The low 
frequency wavenumber 5 signal in m id-latitudes, 
discussed in Section 9, can also be seen here i n the 
tropics in some models such as NCAR. 

The anti- symmetric preci pitation spectra (shown in 
the ATLAS) contain westward propagating non-dis- 
persive modes similar to their symmetric counterparts , 
though with slightly shorter periods in a few models . 
Almost half the models have eastward propagating 
power consistent with the mixed-Rossby-gravity 
(MRG) mode. However, in some of these (ECM-CY32, 
CGAM and UKMO-N48) there is a clear overlapping 
Kelvin mode signature, suggesting that ITCZ preci p- 
itation not entirely symmetric about the equator can 
excite and propagate with Kelvin modes . In fact calcu- 
lations over 3 -years of 6-hourly preci pitation data from 
the NCAR model, which has a double ITCZ, show 
that the time- longitude correlation between precipita- 
tion at the northern and southern hemisphere latitudes 
of maximum zonal average preci pitation is near zero . 
Spectra of outgoing longwave radiation (also shown in 
the ATLAS) are very similar to those of precipitation , 
but lack the background preci pitation spectrum seen in 
a few models. Apparently the more stochastic precipi- 
tation in these models does not prevent their generation 
of equatorial waves. 

Tropical wave spectra are available for Earth climate 
simulations in CMIP3 for several of the APE models, 
documented by Lin et al. (2006). These are GFDL, 
K1JAPAN (named MIROC-medres in CMIP3) , MRI 
and NCAR (run at higher resolution). Of these models, 
only K 1JAPAN has a similar spectrum to observa- 
tions. The others have much weaker power, partic- 
ularly at higher frequencies , and weaker projection 
onto the equatorial wave modes. However, in CMIP3 
each model retains the character of its APE spectrum 
compared to the other models . It is notable that, in 
APE, K ! JAPAN shows perhaps the clearest hierar- 
chical organisation of equatorial convection in the 
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Hovmol 1 er diagrams (Fig. 16). 

It is difficul t to isolate the causes of specific 
behaviour in individual models or groups of models, 
si nee there are multiple differences i n formulation 
between each model. Nor has a clear dependence been 
found on the background flow (mean state). However, 
comparison of ECM-CY 29 and ECM-CY32 indicates 
the effect the parameterizations can have on the equa- 
torial wave spectrum. These two models have the same 
dynam ical core and differ in their parameterization 
suites, including the closure for convective triggeri ng 
and entrainment (Bechtold et al. 2008). The changes in 
ECM -CY32 led to significant improvements in the level 
and realism of forecast tropical variability, documented 
by Bechtold et al., implying that the increased Kelvin 
mode projection in the APE spectrum in ECM -CY32 is 
more "realistic" than ECM-CY29. Wil liamson (2008a) 
studied the dependence of the tropical wave spectrum 
on horizontal resolution and timestep in the NCAR 
model and found that the low resolution (T42) version 
used in APE is far from convergence. In fact the trop- 
ical wave characteristics do not converge for zonal 
wavenumbers less than 16 until T 170 truncation in the 
NCAR model. Compared to T42, the converged spec- 
trum contains relatively more power in lower frequency 
Kelvin modes, rather similar to ECM-CY32 but with 
a stronger background spectrum. However, it is not 
yet known whether a unique converged spectrum will 
emerge in multiple models for the APE CONTROL. 

This analysis has been restricted to parameteriza- 
tion aspects of equatorial variability. Dynamical vari- 
ables would also need to be analysed to fonnally show 
whether the equatorial wave characteristics seen in 
precipitation and OLR in the APE models are indeed 
convectively coupled with the dynamical structures 
predicted by equatorial wave theory. 

1 0.2 Frequency distributions of precipitation 

The frequency distribution of precipitation provides 
infonnation about the extremes simulated by the 
models. In this section we show the frequency distribu- 
tion of 6-hour averaged precipitation in the equatorial 
region, calculated between 1 0°S and 1 0°N to span the 
ITCZ in all the models in Fig. 4. Each column in Fig. 
18 shows a duet of plots. The top row gives the frac- 
tion of time the precipitation is in each 1 mm d- 1 bi n 
ranging from 0 to 120 mm d- 1 • The bottom row gives 
the fraction of time the precipitation is in each 1 0 mm 
d- 1 bin ranging from 0 to 600 mm d 1 • In each plot the 
left-most bin is the fraction of time the precipitation is 
exactly 0, and the right-most bin in the top row is the 
fraction of time the precipitation exceeds 120 mm d 1 • 


Wi 11 iamson (2008a) has shown that the distri butions 
depend on the grid sizes. To el iminate this variation 
the fraction is calcul ated after the model grid data are 
conservativel y mapped to a 5° 1 atitude- longitude grid. 
Plots for the frequency distri butions calculated on the 
original model grids are included in the ATLAS. 

The smallest rates are difficult to discern in Fig. 
18, so Fig. 19 is included to provide details of zero 
and 1 ight precipitation. For this figure the fractions are 
calculated on the origi nal model /data grid, not from the 
data averaged to the 5° grid. For each model this shows 
the fraction of time the 6-hour averaged precipitation at 
grid points between 1 0°S to 1 0°N is zero and the frac- 
tion of time it is positive and less than 0.01 , less than 
O. 1 and less than 1 .0 mm d- 1 *Note that the smallest 
category, 0 < p :S 0.0 1 mm d- 1 , is 1% of the smallest 
bin represented in the first row of Fig. 18. 

The characteristics of trace and no precipitation are 
very different between models. Three models, AGU , 
GFDL and NCAR, rain almost all the ti me, with zero 
preci pitation only 0.005, 0.001 and 0.008% of the time, 
respectively. A few models have zero precipitation for 
a significant fraction of the time: CSIRO, LASG , MR1, 
UKMO-N48 and UKMO-N96 do not rain 27.5, 37.4, 
13.6, 10.8 and 9. 1% of the time , respectively. The 
remaining models range from O. 1 to 4.0%. For the 0 < 
p :S 0.0 1 range most models rain a few percent of the 
time. For the extreme cases, FRCGC rains 68% and 
AG U 0.03% of the time . There is much less model 
spread for 0 < p :S 1.0. Most models are in this range 
around 20% to 30% of the time. The exceptions are 
AG U at 3%, CSIRO at 9% and FRCGC at 80%. In this 
range, however, FRCGC is dominated by the small 0 < 
p :S 0.01 component. 

Figure 18 shows that, away from the very small 
values , the fraction of occurrence decreases mono- 
tonically with increasing rate, except for AGU where 
it peaks at around 7 mm d 1 and CSIRO which has a 
slight peak at around 14 mm day- 1 • ECM-CY29 and 
GFDL show plateaus between 6 and 12 mm d- 1 and 
between 10 and 20 mm d- 1 respectively. 

The models show a wide range of behaviour at large 
precipitation rates. ECM-CY29 has values exceeding 
400 mm d 1 , LASG exceeds 250 mm d- 1 and FRCGC, 
GFDL, K 1JAPAN and MIT reach around 200 mm d- 1 • 
The remaining models have maximum values from 40 
mm d- 1 to 100 mm d- 1 • On the original grids the higher 
resolution models have much higher maximum values, 
with rates greater than 1200 mm d- 1 in FRCGC and 
almost reaching 1200 mm d 1 in ECM-CY 29 (shown 
in the ATLAS). 

As might be expected, the behaviour at large precip- 
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Fig. 18. Fraction of time precipitation is in (top row) 1 mm day- 1 bins ranging from 0 to 120 mm day- 1 , and (bottom 
row) 10 mm day- 1 bins ranging from 0 to 600 mm day- 1 • Grid values between 1 0°S and 1 0°N have been conser- 
vatively averaged to a 5° latitude-longitude fffid . A gap in the curve indicates the fraction is zero for that bin . 


itation rate is consistent with and clearly dominates 
variations in the standard deviation of tropical precipi- 
tation in Fig. 15. The models with a long tail in precip- 
itation frequency and largest maximum values on the 
5° grid are those with a significant background spec- 
trum in Fig. 1 7 and intense quasi- stationary or west- 
ward propagating features in the Hovmciller diagrams. 
The models with a more li mited frequency distribution 
have tropical wave spectra that are restricted to lower 


frequencies and project strongly onto the Kelvin and 
equatorial Rossby modes. Recalling the correlation of 
standard deviation of precipitation with ITCZ width 
noted earlier, a pattern emerges that a narrow single 
ITCZ generally occms i n models that are able to 
generate more intense equatorial precipitation features, 
whereas a broader or double ITCZ generally occurs in 
models with more limited precipitation rates associated 
with domi nant lower frequency variabil ity. 
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Fig. 19. Fraction of time precipitation , p, on the model grid between 10°S and 10°N is (left) equal to zero, and (right) 
in the ranges 0 <p :S 0.01 , 0 < p :S 0.1 and 0 <p :S 1.0 mm day 1 - 


11. Discussion 

The preceding sections describe a range of circu- 
lation statistics from comparison of the APE models 
for the aqua-planet CONTROL SST experiment. The 
conclusions of the individual sections are not repeated 
here and the reader is referred back to each section 
for more detailed discussion . Here, aspects of the 
aqua-planet simulations worthy of particular note are 
discussed and more general conclusions are drawn 
from the model comparison . 

A multi-model mean and standard deviation have 
been computed , as a means of showing variation 
between the models in a concise form. The multi- 
model mean is not intended to be a reference solution 
for the CONTROL aqua-planet, whose circulation 
cannot be precisely known. The aqua-planet climate 
is very earth-like in structure, but the SST profile , 
its zonal symmetry, the absence of orography and 
perpetual equinoctial insolation combine to create a 
stronger circulation, with both a stronger westerly jet 
and tropical meridional circulation . 

In comparing atmospheric model simulations on an 
aqua-planet , it is not clear a priori whether the range of 
model behaviour will be larger or smal ler than for an 
equivalent comparison of simulations of Earth climate. 
On the one hand, removing land, orography and zonal 
SST variations results in a simpler circulation, which 
might lead to an expectation of reduced model spread. 
Conversely, the lack of zonally asymmetric forcing 
(additional to the diurnal cycle of insolation) impl ies 
that all longitudinal structure and variability must 


be internally generated by the model dynamics and 
parameterizations and by interactions between them. 
This might lead to an expectation of increased model 
spread and was an important motivation in the design 
of the idealised aqua-planet configuration by Neale and 
Hoskins (2000a). 

Indeed , the comparison does reveal a wide range of 
model behaviour in most of the atmospheric circula- 
tion statistics considered . However, certain parameters 
are more constrained , with relatively small differ- 
ences between the models. Generally, it is aspects 
of the circulation determined mainly by the resolved 
dynamics that display less variation , whereas aspects 
that involve interaction between dynamical and param- 
eterized moist processes appear to be most variable 
between the models. To illustrate this for the most 
basic statistic, namely the time average zonal average 
tropospheric state, the balanced component (tempera- 
ture and zonal wind) varies relatively little between the 
models, while the tropical divergent circulation varies 
significantly, associated with the strength, width and 
vertical structure of the equatorial ITCZ. The trade 
wind inversion and longwave cooling in the sub-tropics 
associated with boundary layer cloud also vary signifi- 
cantly between the models. Furthermore , even for the 
dynamical circulation, quantities influenced by budget 
constraints have smaller variation . For example, the 
transient eddy covariances (poleward fluxes of heat and 
momentum) show significantly less vaiiation between 
the models than the transient eddy variances (including 
kinetic energy), which are particularly sensitive to 
model resolution and numerical damping. 


52 


Journal of the Meteorological Society of Japan 


Vol. 91A 


In contrast to the tropospheric mean state, the lower 
stratosphere varies greatly between the models, due, at 
least i n part, to the proximity and differing numerical 
treatments of the upper boundary. The historical uncer- 
tainty in both polar and tropical lower stratospheric 
temperatures persists in APE, with a larger range of 
tropical tropopause temperature and vert ical structure 
than in a recent AMIP comparison (part of CMIP3). 
This may reflect the i nclusion in APE of a number of 
development models, compared to the more estab- 
lished models in AMIP and other intercomparisons. 

All aspects of tropical precipitation vary greatly 
between the models. In the time average, the ITCZ 
appears as either a single peak on the equator or a 
double peak with a relative minimum at the equator 
in individual models, i n common with si mulations of 
Earth climate (Lin 2007, and references therein). This 
behaviour is sensitive to the latitudi nal profile of SST, 
with evidence of a transition from a single equato- 
rial ITCZ to a widel y separated double ITCZ as the 
tropical SST profile flattens, as discussed by Neale 
and Hoskins (2000b) and Williamson et al. (2013). 
In single model sensitivity studies (using the NCAR 
model), Wi lliamson and Olson (2003) and Williamson 
(2008a) have shown that the ITCZ structure is sensi- 
tive to the convective parameterization, to detai Is of its 
implementation and to its interaction with the model 
dynamics . 

The spatio-temporal variability of tropical preci pi- 
tation displays varying degrees of projection onto the 
theoretical equatorial wave modes . In most models 
with a strong projection, convective features propagate 
predominantly eastward, but in a few models westward 
propagation dominates . Eastward propagation gener- 
ally projects onto the non-dispersive Kelvin modes 
and has remarkably constant phase speed among the 
models. Westward propagation is less dispersive than 
the theoretical equatorial Rossby modes . Several 
models display the hierarchical convective organisa- 
tion seen in Earth observations (N akazawa 1 988), with 
smaller scale westward propagation appearing within 
an envelope of eastward propagation . This is particu- 
larly true of the higher resolution models that have grid 
spacings of order 100 km. 

A relationship exists between the tropical mean state 
and variability among the APE models, consistent with 
the results of previous studies. A nan-ow single ITCZ 
tends to occur in models with more intense equatorial 
precipitation features, and a wider or double ITCZ i n 
models dominated by lower frequency tropical waves 
and weaker precipitation variability. The reason for 
this relationship is unclear, beyond a scaling of precip- 


itation variabi lity by the equatorial zonal average, or 
does the relationship appear predictive of the transition 
to a widely separated double ITCZ when the tropical 
SST profile is flattened in the QOBS and FLAT SST 
experi ments, discussed by Wi 1 liamson et al. (201 3). 
Two of the models with a narrow single ITCZ in the 
CONTROL experiment are among those most strongly 
split in QOBS. 

The global energy balance is smprisi ngly uncon- 
strained between the APE models , because it is one 
parameter that is strongly constrained in si mulations 
of Earth cli mate. The global average TOA net flux is 
of course unknown for the aqua-planet, but it has a 
range of 30 W m 2 in the APE models, dominated by 
variations in shortwave reflection by clouds. In most of 
the models the TOA net flux is relatively insensitive to 
the aqua-planet SST profile, so the main contributing 
factors to this large range appear to be a lack of opti- 
misation for present-day Earth cl imate in a majority 
of the models and, more speculatively, the presence of 
differing biases in TOA net flux over land in the models 
that are optimised. To resolve this uncertainty, it would 
be necessary to quantify the global energy balance both 
in aqua-planet and AMIP simulations of Earth climate 
for individual models. This illustrates the intention of 
Neale and Hoskins (2000a) that APE be considered 
one component of a modelling hierarchy. The detailed 
aqua-planet circulation is ultimately unknown, though 
its overall character is predicted by theory and ideal- 
ised models. Aqua-planet simulations are therefore 
most usefully compared and assessed in relation to 
simulation of the known Earth climate. 

The large differences seen in global average TOA net 
energy flux among the APE models also extend to its 
latitudinal distribution, forcing large differences, with 
a standard deviation of around 1 PW in mid-latitudes , 
in the total poleward energy transport of the atmo- 
sphere-ocean system required for balance. This implies 
that significantly different equilibrated climates would 
result if the APE models were coupled to either ther- 
modynamic slab oceans or dynamic ocean models . 

It is commonly thought that much of the difference 
between atmospheric GCM simulations is due to the 
sub-grid scale parameterizations. The results of APE 
are not sufficiently generic to allow us to comment 
on systematic impacts of particular modelling choices 
such as parameterization scheme or physics-dynamics 
coupling across the APE models. This requires experi- 
ments specifically targeted to such questions in a more 
constrained model ling environment. However, part of 
the difference in APE might be primarily due to the 
truncation error of the numerical approximations to the 
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resolved fluid flow component of the models, especially 
as many of the model resolutions are rather coarse, i n 
the area of T42 spectral truncation or 2.5° grids. These 
do not ful ly capture barocli nic wave development and 
differ significantly from the same models at higher 
resolution (Jablonowski and Will iamson 2006). This 
will influence the eddy variances and covariances in 
particular. However, even in this case, parameteriza- 
tion forcing may play a role. 

Resol ution will also influence phenomena that arise 
fundamentally from interaction between the model 
dynam ics and physical parameterizations, particularly 
if the response of either model component to changes 
in the other is non -linear. In the APE si mulations this 
most clearly appl ies to the convectivel y coupled equa- 
torial waves. Williamson (2008a) has shown that in the 
NCAR CAM the tropical wave characteristics depend 
on resol ution and for zonal wavenumbers less than 16 
they do not converge unti 1 T 170 truncation, while T42 
is common in the APE experiments here. This lack 
of convergence at the resolutions that APE models 
have been run contributes to the difference among the 
models. 

It is important, therefore, that convergence of the 
APE mean state and variability with increasing reso- 
1 ution is studied in a number of models. It is conceiv- 
able that individual models could converge towards a 
unique climate but, equally, they might converge to 
significantly different climates, because of their wide 
range of parameterizations and, possibly to a lesser 
extent, their dynamical approximations . 

A future possibility might be the application of 
global cloud resolving models (CRMs) to the aqua- 
planet configuration , as computational capacity 
continues to increase. CRMs do retain parameteriza- 
tions, of radiation , turbulence and cloud microphysics, 
but the more empirically based closure assumptions 
of convective parameterizations in large scale models 
are not required . If global CRMs were to converge to 
a single climate, more confidence would be gained 
that this was indeed the aqua-planet climate . However, 
confidence would also be required in the CRMs them- 
selves, for example by accurate simulation of Earth's 
global climate or at least the particular phenomena 
relevant to an aqua-planet. Progress has already been 
made in this regard by coarser resolution cloud-system 
resolving models (CSRMs, one of which is included 
in APE) . In a number of local case studies organised 
within the GEWEX Cloud Systems Studies (GCSS), 
CSRMs have proved successful in simulating observed 
convection more accuratel y than have parameteriza- 
tions extracted from large scale models . 
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